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INTRODUCTION 


The  formalities  associated  with  this  contract  were  completed  during  October  1986, 
and  the  effective  starting  date  of  the  project  is  October  1,  1986.  The  present  report  covers 
the  period  from  the  starting  date  until  October  1,  1987. 

The  main  objective  of  this  project  is  to  determine  conditions  of  exposure  of  human 
lymphocytes  in  vitro  to  continuous  wave  (CW)  or  pulsed  wave  (PW)  2450  MHz  microwaves, 
which: 

•  affect  the  course  of  spontaneous  or  phytohemagglutin  in  (PHA) —  induced  lympho- 
blastoid  transformation,  and 

•  may  induce  numerical  and/or  structural  aberrations  in  chromosomes  of  cells  undergo¬ 
ing  division  following  transformation. 

In  order  to  provide  reliable  quantitative  data  on  exposure,  a  system  with  automated 
dosimetry  was  developed,  and  tested  for  biocompatibility  with  lymphocyte  cultures  in 
vitro  for  periods  of  up  to  120  hrs.  (5  days).  A  method  for  quantitation  of  lymphocyte 
transformation  based  on  image  analysis  was  developed  and  tested.  Experiments  on  effects 
of  exposure  on  spontaneous  transformation  were  carried  out.  For  convenience  and  clarity 
the  work  performed  is  described  under  two  separate  headings:  (1)  exposure  system  and 
dosimetrv,  and  (2)  studies  on  lymphocyte  cultures. 

f 

EXPOSURE  SYSTEM  AND  DOSIMETRY 
Description  of  the  Exposure  System 

The  experiments  planned  in  this  project  necessitated  the  design  and  assembly  of  an 
exposure  system,  which  would  meet  several  engineering,  dosimetric  and  biological  require¬ 
ments.  The  system  should  provide: 

1)  an  environment  compatible  with  lymphocyte  (tissue)  culture  in  vitro  in  respect  to 
relative  humidity  (about  100%)  CO2  content  in  air  (5%)  and  temperature  (37°C);  an 
additional  requirement  for  temperature  regulation  consisted  in  providing  for  heating 
and  equilibrating  temperature  of  the  samples  up  to  39°C  by  microwave  or  conventional 
heating; 
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2)  the  possibility  of  exposing  samples  (preferably  exposure  of  multiple  samples  under 
identical  condition)  to  CW  or  PW  2450  MHz  with  a  high  ratio  of  peak  to  average 
power,  and  the  possibility  of  simultaneous  sham  exposure  of  control  samples  under 
strictly  comparable  conditions,  where  the  only  difference  between  exposed  and  sham 
exposed  samples  would  consist  in  the  presence  (or  absence)  of  microwave  absorption 
in  the  sample; 

3)  continuous  measurements  and  registration  of  the  temperature  of  exposed  and  sham 
exposed  samples  during  exposure  to  provide  a  temperature  profile  over  time  history 
of  the  samples,  and  data  for  computation  of  the  specific  absorption  rate  (SAR)  and 
specific  absorption  (SA)  for  each  of  the  experimental  conditions  tested; 

4)  continuous  monitoring  of  the  wave  form  and  repetition  rate  in  the  case  of  PW  expo¬ 
sures. 

The  first  and  second  requirements  were  met  by  placing  two  shorted  rectangular  S-band 
waveguides  300  mm  long  (Omega  Laboratories  Models  6101)  in  a  conventional  CO2  tissue 
culture  incubator  (Forma  Scientific)  thermostated  at  37°C.  A  styrofoam  block  placed  in  the 
waveguide  serves  as  support  for  the  sample  holder,  and  ensures  that  the  sample  is  centrally 
placed  in  the  waveguide  in  the  same  location  in  successive  experiments.  A  small  hole 
drilled  in  the  waveguide  allows  the  introduction  of  a  non-perturbing  temperature  probe. 
The  hole  is  drilled  near  the  guide  wall,  at  the  center  of  the  shorter  dimension,  this  is  a 
low  electric  field  point  and  the  hole  produces  minimal  disturbance.  One  of  the  waveguides 
serves  for  sham  exposures,  the  other  one  is  connected  through  a  matched  coaxial  cable  to 
the  microwave  power  source.  This  consists  of  several  elements.  The  output  from  a  CW 
Hewlett-Packard  Model  8616A  oscillator  feeds  a  Hughes  Model  1177H  travelling  wave  tube 
amplifier.  Pulsed  or  amplitude  modulated  waveforms  can  be  obtained  through  the  use  of 
a  Hewlett-Packard  Model  8403A  p-i-n  modulator,  for  high  peak  power  pulsed  exposures, 
the  signal  can  be  further  amplified  with  a  1  kW  amplifier  (MCL,  La  Grange,  IL).  The 
waveguide  receives  its  power  through  an  isolator  and  a  dual  directional  coupler  that  allows 
forward  and  reflected  powers  to  be  measured  with  Hewlett-Packard  Model  432A  power 
meters  and  waveforms  to  be  monitored  using  a  Phillips  PM  3240  oscilloscope.  The  amount 
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of  power  supplied  to  the  waveguide  is  controlled  by  the  output  from  the  oscillator  and  the 
voltage  of  the  1  kW  amplifier,  and  monitored  by  measurements  of  the  forward  and  reflected 
powers,  and  by  the  readings  on  the  oscilloscope.  The  main  function  of  the  oscilloscope  is, 
however,  to  monitor  waveforms.  These  means  of  control  and  monitoring  are  only  auxiliary, 
as  the  conditions  of  exposure  are  characterized  by  dosimetry  based  on  measurements  of 
temperature  and  of  SAR  within  the  sample,  as  specified  by  the  3rd  requirement  the  system 
should  meet. 

The  dosimetric  part  of  the  system  consists  of  two  non-perturbing  temperature  probes, 
connected  to  a  Hewlett-Packard  Model  59306A  relay  activator,  which  is  itself  connected 
to  a  Keithley  Model  192  digital  voltmeter  (DVM).  Both  the  relay  activator  and  the  DVM 
are  under  control  of  a  Hewlett-Packard  Model  86  desk-top  computer  through  a  HP1B 
(IEEE-488  bus).  Initially  a  Narda  Model  801  IB  non-perturbing  double  temperature  probe 
was  used,  later  two  Vitek  Model  101  probes  were  substituted  for  the  Narda  probe.  The 
Narda  probe  requires  cumbersome  calibration  against  a  NBS  standard  and  introduction  of 
correction  factors  to  relate  the  current  to  temperature,  a  linear  relationship  being  main¬ 
tained  over  limited  temperature  ranges.  Vitek  probes  are  much  easier  to  calibrate,  are 
more  stable,  and  the  current-temperature  relationship  is  linear  over  the  temperature  range 
of  interest. 

Under  computer  control  the  temperatures  of  exposed  and  sham  exposed  samples  are 
recorded  sequentially  and  stored  in  memory  (records  are  maintained  by  keeping  a  library 
of  micro-floppy  disks).  Since  the  temperature  of  the  sham  exposed  sample  is  essentially 
constant  in  function  of  the  stability  of  the  temperature,  of  the  conventional  incubator  it 
can  be  monitored  less  frequently  than  the  temperature  of  the  exposed  samples.  For  conven¬ 
tional  thermal  exposures  the  temperature  can  be  elevated  by  increasing  the  temperature 
of  the  incubator,  maintaining  control  samples  in  another  CO2  incubator.  For  short  term 
conventional  heating  experiments  a  Precision  (GCA  Corporation)  water  bath  can  be  used 
inside  the  incubator. 

The  sample  holder  has  to  be  a  tissue  culture  vessel  made  of  materials  non-toxic  to 
cells  in  vitro  and  transparent  to  2450  MHz  microwaves.  The  geometry  of  the  container  has 
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to  allow  for  uniform  energy  deposition  within  the  sample.  These  considerations  limited  the 
choice  of  possible  sample  holders.  Previous  experiences  with  tissue  culture  exposures  in 
a  S-band  waveguide  demonstrated  that  the  use  of  T~flasks  or  Petri  dishes  leads  to  signifi¬ 
cant  nonuniformities  in  energy  deposition.  Two  types  of  sample  holders  were  tested,  both 
rectangular  tissue  culture  Lab-Tek  (Miles  Laboratories)  chamber  slides,  a  two  chamber 
model  N04S02,  and  a  four  chamber  model  No4804.  Analysis  of  SAR  distribution  based 
on  measurements  of  temperature  (see  below)  in  slide  chambers  containing  tissue  culture 
medium  demonstrated  that  the  energy  distribution  within  and  between  chambers  of  the 
two  chamber  model  is  nonuniform,  the  difference  between  chambers  being  about  50%.  In 
the  four  chamber  model  measurements  demonstrated  that  differences  in  energy  deposition 
remain  within  10%  of  an  average  value.  Therefore,  the  four  chamber  model  was  selected 
for  use  in  experiments.  Three  of  the  chambers  contain  tissue  culture  medium  and  cells, 
the  fourth  contains  medium  only.  A  hole  is  drilled  in  the  cover  of  the  chamber  and  a  tem¬ 
perature  probe  is  introduced.  In  this  way  each  experiment  can  be  carried  out  in  triplicate 
on  simultaneously  exposed  samples  with  simultaneous  continuous  dosimetric  control. 

Dosimetry  and  Temperature-Time  Anaysis 

During  an  experiment  the  temperature  in  the  exposed  and  the  sham-exposed  samples 
is  recorded  at  regular  intervals,  with  a  minimum  interval  of  less  than  1  s.  The  temperature 
is  recorded  before  the  beginning  of,  during  and  after  the  exposure.  The  “on’1  and  “off* 
times  of  exposure  are  recorded  on  the  computer  by  the  operator  using  a  soft-key  interrupt 
capability.  At  the  conclusion  of  a  run  the  T(t),  i.e.,  temperature  (T)  versus  time  (f) 
behavior  is  analyzed  to  determine  SAR. 

If  the  SAR  is  5  (\Vkg~1)  and  the  specific  heat  of  the  sample  is  C(Jkg“1K“1),  then  the 
rate  of  heating  during  microwave  exposure,  for  small  temperature  variations  over  which  C 
can  be  assumed  to  be  constant,  is 


If,  prior  to  deliberate  heating,  a  sample  is  below  its  equilibrium  temperature,  Ttq ,  with  its 
surroundings,  then  for  small  temperature  differences  from  equilibrium,  the  rate  of  temper- 
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ature  change  is  of  the  form 

(f)+="<r--T)'  <2> 

Similarly,  if  after  an  exposure  is  terminated  a  sample  cools  the  rate  of  temperature  change 
will  be 

<3> 

Eqs.  (2)  and  (3)  are  equivalent  statements  of  Newton’s  law  of  cooling.  In  general,  when 
microwave  power  is  applied  to  a  sample  that  was  not  in  thermal  equilibrium  with  its 
surroundings  at  the  start  of  exposure: 

dt  c+a^T'q  ^  ^ 

The  solution  to  Eq.  (4)  can  be  written  in  the  form 

T  -  To  =  +  Tcq  -  To)  (1  -  e~a<)  (5) 

where  T0  is  the  control  temperature  at  t  =  0.  Either  Eq.  (4)  or  Eq.  (5)  can  be  used  to 
determine  the  SAR.  From  Eq.  (4), 


so  that  determination  of  the  heating  rate  before  and  during  application  of  microwave  power, 
or  during  and  after  the  application  of  microwave  power,  determines  the  SAR. 

If  the  sample  is  in  thermal  equilibrium  before  exposure  begins,  Eq.  (5)  reduces  to  the 
simpler  form 

T~To  =  ^(1~e0,)  (7) 

which  can  be  used  to  find  S  from  T(t).  In  general,  the  approach  using  Eq.  (6)  is  preferable, 
since  for  times  near  t  =  0,  the  exponential  can  be  expanded  in  quadratic  form. 
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In  practice,  this  analysis  is  performed  automatically.  At  the  end  of  a  run,  the  point  of 
inflection,  or  turning  point,  in  T(t)  is  found  numerically.  The  temperature  profile  to  the 
left  and  right  of  this  point  is  fitted  by  least-squares  to  a  linear  or  quadratic  function  and 
the  change  in  slope  at  the  beginning  of  the  exposure  period  yields  the  SAR. 

In  practice  temperatures  in  the  exposed  sample  are  recorded  each  1,  2  or  3  seconds, 
the  readings  from  the  sham-exposed  sample  being  taken  at  10  times  longer  intervals.  In 
experiments  in  which  the  sample  is  microwave  heated  for  longer  periods  at  a  predetermined 
elevated  temperature,  measurements  are  made  less  frequently  once  equilibrium  is  reached, 
and  serve  only  as  control  of  the  stability  of  exposure  conditions. 

Measurements  made  in  air  or  water  in  the  incubator  outside  the  waveguides,  demon¬ 
strated  oscillations  of  ±0.2°C  in  dependence  from  the  on-off  cycle  of  heating  controlled 
by  the  thermostat.  Over  a  24  hour  period  such  oscillations  may  reach  ±0.5°C.  Measure¬ 
ments  in  medium  in  chamber-slides  placed  within  the  waveguides  demonstrated  that  these 
oscillations  are  attenuated,  the  metal  structure  of  the  waveguide  acting  as  a  heat  sink  or 
source .  Over  a  24  hour  period  the  oscillations  do  not  exceed  0.05°C,  and  usually  remain 
within  0.02°C.  The  sensitivity  of  the  temperature  probes  is  nominally  0.01°C,  however 
comparison  of  successive  readings  permits  extrapolation  to  0.005°C.  Based  on  data  in  the 
literature  the  biological  endpoints  examined  are  not  affected  by  temperature  excursions 
within  0.5°C. 

The  software  written  for  reading  temperatures  from  two  Narda  probes  and  for  “smooth¬ 
ing  and  plotting  of  experimental  data",  i.e.,  for  SAR  computation  is  enclosed  as  Appendix 
1.  The  software  has  provision  for  temperature  readings  to  be  taken  by  Narda  or  Vitek 
probes. 

The  dosimetric  system  described  above  can  be  used  in  experiments  on  the  effects  of 
heating  not  only  by  microwaves  but  also  by  any  means  (eg.  conventional,  ultrasound). 
Because  of  its  wide  applicability  the  system  is  the  subject  of  two  brief  communications: 
one  at  the  IRPA-7  Congress  in  April,  1988  in  Sydney  Australia  (Appendix  2),  and  a  second 
one  at  the  10th  Annual  Meeting  of  the  BEMS  in  June,  1988  in  Stamford,  CT  (Appendix  3). 
Finally,  a  detailed  description  has  been  submitted  and  accepted  for  publication  in  Health 
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Physics  (Appendix  4). 


STUDIES  OF  LYMPHOCYTE  CULTURES 
Collection  of  Blood  Samples  and  Separation  of  Lymphocytes 

Approval  for  the  use  of  human  blood  was  obtained  November  19,  1986  from  the  Human 
Subjects  Research  Review  Board,  Office  of  the  Surgeon  General,  Department  of  the  Army. 
Volunteers  are  informed  about  the  nature  of  the  study,  potential  risks  and  the  potential 
benefit  derived  from  the  study,  and  sign  the  informed  consent  form,  as  approved  by  the 
HSRRR.  20ml  of  blood  is  withdrawn  by  venipuncture  under  aseptic  conditions  into  two 
heparinized  (143  USP  units)  sterile  evacuated  blood  collection  tubes  (Vacutainer  No.  6480, 
Becton  and  Dickinson).  The  blood  is  transported  to  the  laboratory  at  room  temperature. 
Within  one  hour  from  withdrawal  mononuclear  cells  (lymphocytes  and  monocytes)  are 
separated  under  sterile  conditions  from  the  whole  blood  by  gradient  centrifugation  in 
Lymphocyte  Separation  Medium  (Boehringer,  Mannheim,  West  Germany)  according  to 
the  manufacturers  instructions.  Following  separation  the  cell  suspension  is  w’ashed  twice 
in  RPMI  1640  cell  culture  medium  (Gibco).  An  aliquot  is  withdrawn  for  the  determination 
of  the  cell  concentration,  cell  viability  and  cytological  (cytocentrifuge)  preparations,  the 
remainder  is  used  for  establishing  cultures. 

The  cells  are  counted  in  a  bright  line  hemocytometer  (Neubauer  type,  Spencer).  Vi¬ 
ability  is  tested  by  the  Trypan  blue  exclusion  test,  neutral  red  stain  for  the  vacuole,  and 
Janus  green  stain  for  mitochondria,  Cytocentrifuge  (Cytospin  2  cytocentrifuge,  Shandon) 
preparations  are  air  dried,  fixed  in  methanol,  stained  by  a  combined  Wright-Giemsa  stain 
and  cell  morphology  is  analyzed  under  a  microscope  (see  below). 

The  usual  yield  of  cells  ranges  from  1.107  to  3xl0‘.  The  suspension  consists  of 
90-95%  lymphocytes,  the  remainder  of  cells  being  monocytes  and  occasional  (below  1%) 
granulocytes.  The  viability  of  cells  ranges  from  96.5  to  99.5%,  usually  is  about  98 

Cell  Cultures 

Four-chamber  culture  chamber  slides  are  filled  under  sterile  conditions  with  1  ml  of 
chromosome  medium  1A  (Gibco)  with  or  without  phytohemogglutin  (PHA)  depending 
on  the  protocol  appropriate  for  the  particular  experiment.  106  cells  axe  added  to  3  of  the 

7 


chambers,  the  4th  contains  medium  only.  A  0.2  mm  hole  is  drilled  over  the  chamber,  which 
does  not  contain  cells,  for  the  introduction  of  the  temperature  probe.  For  experiments  on 
spontaneous  transformation  a  set  of  experimental  samples  consists  of  four  slides:  one  with 
PHA  (5  or  10  mg/ml)  positive  control  incubated  in  a  conventional  incubator,  and  three 
without  PHA.  Out  of  these  one  is  incubated  in  a  conventional  incubator  at  37°C,  one  is 
microwave-exposed  and  one  is  sham  exposed. 

Cell  Harvest 

Following  incubation  the  cultures  without  PHA  are  mixed  within  the  culture  chamber 
using  a  Pasteur  pipette,  and  transferred  to  microfuge  tubes.  The  cell  count  is  established, 
viability  is  tested  and  cvtocentrifuge  preparations  are  made  as  described  above. 

In  the  case  of  cultures  incubated  with  PHA  10/il  of  colcemide  solution  (10  mg/ml, 
Gibco)  are  added  to  each  chamber  2  hours  before  harvest.  At  harvest  the  cultures  are 
transferred  to  microfuge  tubes,  spun  down  at  103  rpm  for  8  min  in  an  Eppendorf  micro 
centrifuge,  and  the  supernatant  is  discarded.  1  ml  1%  sodium  citrate  solution  prewarmed 
to  37°C  is  added  with  vortexing.  The  cells  are  kept  in  this  solution  for  8  min,  spun  down 
under  the  same  conditions  as  above,  and  the  supernatant  is  discarded.  Fixative  (cold 
glacial  acetic  acid:  methanol,  1:3  vol/vol)  is  added  dropwise  with  vortexing.  After  at  least 
20  min  the  fixative  is  changed  once  or  twice.  A  droplet  of  the  cell  solution  in  fixative  is 
placed  on  a  microscope  slide  wet  with  20%  ethanol,  and  the  slides  are  flamed.  Following 
this  the  slides  are  air  dried,  and  stained  with  4solution  in  phosphate  buffer  pH  6.8  for  8 
min  or  treated  with  trypsin  and  stained  with  Giemsa  according  to  the  seabright  procedure 
(G-banding). 

Analysis  of  Microscope  Preparations 

Cytocentrifuge  slides  are  examined  under  a  lOOx  oil  immersion  objective  and  the 
cell  and  nucleus  areas  are  determined  using  an  “Optomax”  image  analyzer.  The  original 
software  supplied  by  the  manufacturer  was  rewritten.  The  slides  are  coded  and  the  person 
who  makes  the  measurements  is  not  aware  of  the  treatment  (exposed,  sham-exposed  or 
conventional  incubator).  The  code  number  of  the  slide,  the  area  of  the  cell,  the  area  of  the 
nucleus,  the  nucleus-cytoplasmic  ratio  and  a  symbol  for  morphological  features  of  the  cell 
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(S  for  small  lymphocyte,  1  for  an  intermediate  form  and  B  for  a  lymphoblastoid  cell)  are 
recorded  and  stored  on  floppy  disk. 

Chromosome  preparations  are  examined  under  a  100 X  oil  immersion  objective,  and 
the  number  and  structure  of  chromosomes  is  noted.  In  this  part  of  the  study  chromosome 
preparations  were  made  only  from  positive  control  (with  PHA)  samples,  no  studies  on 
microwave  exposed  samples  were  carried  out,  as  the  number  of  cell  divisions  following 
spontaneous  transformation  is  too  low  for  any  meaningful  analysis. 

RESULTS 

A  series  of  preliminary  experiments  with  samples  incubated  with  PHA  for  48  or  72 
h  demonstrated  that  the  tissue  culture  technique  and  harvesting  technique  as  modified 
for  4-chamber  slides  is  adequate.  Lymphoblastoid  transformation  of  85-95%  of  cells  was 
obtained  at  PHA  concentration  of  10  mg/ml  in  medium  after  72  h  of  incubation  (manu¬ 
facturer  recommended  concentration).  At  suboptimal  concentrations  of  2.5  or  5  mg/ml. 
transformation  was  obtained  in  about  20  or  40-50%  respectively.  At  5  mg/ml  concentration 
the  mitotic  index  is  about  10%.  In  view  of  these  results  it  was  concluded  that  the  culture 
technique  and  conditions  in  the  exposure  system  axe  compatible  with  the  requirements  of 
human  lymphocytes  in  vitro.  Cell  viability  following  incubation  remained  within  the  same 
range  as  in  the  initial  blood  sample. 

Experiments  with  spontaneous  transformation  (incubation  without  PHA)  were  carried 
out  at  nonheating  and  heating  (0.5°,  1.0°  or  2°C)  levels  of  exposure  both  with  continuous 
waves  and  pulsed  waves  (P\V).  At  heating  levels  of  a  2°C  increase  after  5  or  4  day  exposure 
a  significant  decrease  of  cell  viability  occurs.  Cell  debris  was  noted  in  preparation.  Thus 
experiments  with  heating  at  39°C  or  more  were  discontinued  as  cell  death  introduces 
confounding  factors. 

Following  C W  exposures  at  nonheating  levels  or  heating  by  0.5°C  no  differences  be¬ 
tween  sham-exposed  samples  or  control  samples  incubated  in  a  conventional  incubator 
were  seen.  Heating  by  1.0°C  has  a  slight  effect  in  enhancing  spontaneous  transformation. 

Exposures  to  pulsed  2450  MHz  at  both  nonheating  and  heating  levels  enhances  spon¬ 
taneous  transformation  when  compared  to  sham-exposed  samples  or  samples  incubated  in 


9 


a  conventional  incubator. 


The  results  mentioned  above  are  preliminary.  The  statistical  significance  of  the  results 
is  now  being  analyzed  using  a  multivariate  analysis  program  (Statgraph).  Spontaneous 
transformation  varies  individually  with  donors.  Therefore  the  sample  incubated  in  a  con¬ 
ventional  incubator  provides  a  baseline  with  which  transformation  in  the  sham-exposed  or 
exposed  sample  are  compared.  The  need  for  comparison  of  results  from  each  individual 
set  of  samples  from  the  same  donor  with  a  set  from  another  makes  the  statistical  analysis 
more  complex,  as  the  significance  of  the  magnitude  of  differences  from  the  base-line  in 
each  set,  depending  on  exposure  conditions  has  to  be  established.  It  can  be  concluded, 
however,  that  measurements  of  the  cell  area  provide  a  good  index  of  the  dynamics  of  the 
transformation  process. 
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program  for  collecting  data 
OPTION  BASE  ! 

DISP  "ARE  YOU  USING  VITEK  OR  NARDA  PROBES, TYPE  V  OR  N" 

INPUT  PRBTYPS 

PR0BCI=9.6184  9  PR0BC2=41  .459  !  Calibration  factors  for  Naraa  proDe  1 
PRGBC3  =  S .  91 21  £  PR0BC4=42 . 38b  !  Calibration  factors  for  Narda  prone  2 
IF  PRBTYP$="N”  THEN  GOTO  90  !  Probe  type  Narda  -  don't  reset  caiib.  factors 
PR0BC1  =  1 00  £  PRQBC2  =  0  !  Calibration  factors  for  Vitek  probe 
PR0BC3=  1 00  PR0BC4  =  0  1  Calibration  factors  for  Vitek  probe 
ON  KEY*  1  , "POWER  ON"  GOSUB  1380 
)  ON  KEY*  4, "NEW  OBS  INT"  GOTO  1100 
J  ON  KEY*  7. "POWER  OFF"  GOSUB  1390 
)  MASS  STORAGE  IS  ":D500" 

)  DIM  ATEMP( 1000) , ATIME < 1 000 > , BTEMP( 1 000  > . BT IME( 1 000 ) , PONT ( 50 ) ,PGFFT<50) 

)  DIM  ANUM(20> ,BNUM(20> ,IPON(20> ,1P0FF<20> 

)  CLEAR 

)  !  BEEP  40,100  »  WAIT  50  #  BEEP  40,100  $>  WAIT  50  *  BEEP  40,100  «  WAIT  50  I*  BE 
80,150 

)  DISP  "ENTER  DAY .MONTH, YEAR" 

)  INPUT  EXPDATES 
)  EXPDATE=0 

)  DISP  "ENTER  SAMPLE  IDENTIFIER" 
j  INPUT  SAMPLES 
)  EXPDATAS=EXPDATES4SAMPLES 

)  DISP  "ENTER  TIME  OF  DAY  -  24  HOUR  CLOCK"  <s  INPUT  TODS 
)  ON  ERROR  GOTO  300 
)  CREATE  EXPDATAS6" : D50 1 ",  1 .32000 
)  ASSIGN*  1  TO  EXPDATAS&" : D50 1 " 

)  ’  PRINT*  1  ;  EXPDATES 
)  1  PRINT*  1  ;  SAMPLES 
j  GOTO  390 

)  IF  ERRN  =128  THEN  DISP  "DISC  IS  FULL  -  INSERT  A  NEW  DISC  AND  START  AGAIN" 

J  IF  EKKN  =128  THEN  GOTO  3h0 
)  IF  ERRN  =83  THEN  GOTO  350 

)  DISP  "AM  ERROR  HAS  OCCURRED  -  PRESS  CONTINUE  TO  START  AGAIN" 

)  PAUSE  8  GOTO  150 

l  DISP  "THIS  DA i A  FILE  NAME  ALREADY  EXISTS  -  ENTER  ANOTHER  SAMPLE  NUMBER" 

J  INPUT  SAMPLES 
)  EXPDATAS=EXPDATES&SAMPLE$ 

)  GOTO  250 

!  KK  =  0  ?  Counter  for  new  observation  times  for  probes  A  and  B. 

J  BlN=0  !  Holds  previous  number  of  observations  from  probe  A. 

I  NDP I S  =  0  !  Counter  for  number  of  data  points  recorded. 

)  CLEAR 

l  DISP  This  program  uiili  store  1  000  temps,  and  times  for  ootn  proDes." 

)  DISP  "ENTER  TOTAL  OBSERVATION  TIME  IN  SECONDS." 
i  INPUT^OBSTIME 

1  DISP  "Probe  A  should  be  in  the  exposeo/exper imental  culture." 

I  DISP 

DISP  "ENTER  OBSERVATION  INTERVAL  FOR  PROBE  A  IN  SECONDS." 
i  INPUT  NA 

IF  NA>OBSTIME  THEN  520 
i  GOTO  540 

,  rriTn  se^ected  ods  .  interval  >  the  total  obs.  time;  try  again." 

1  bU  I  U  4  J  U 

AjQBS=INT  (OBSTIME/NA) 

NA=NA-1000 

hold=ajgbs+bin 

IF  H0lD> 1 00U  THEN  590 
GOTO  610 

^°U  ^ave  requested  too  many  readings  of  probe  A;  try  again." 

b  J  I  u  A  J  Li 


610  DISP  "The  observation  interval  for  probe  B  should  be  an  integer  multiple" 

S20  DISP  "of,  and  ideally  <  or  =  10  times  that  for  probe  A." 

630  DISP 

640  DISP  "ENTER  OBSERVATION  INTERVAL  FOR  PROBE  B  IN  SECONDS." 

650  INPUT  NB 

660  BJ0BS= INT  (OBSTIME/NB) 

670  NB=NB*1 000 

630  IF  FP  <NB/NA>*0  THEN  700 
690  GOTO  720 

700  DISP  "The  obs.  time  for  B  is  not  a  multiple  of  that  for  A;  try  again." 

710  GOTO  640 

720  REMOTE  502  !  DVM  to  remote. 

730  CLEAR  9  KEY  LABEL 

740  REMOTE  507  !  Rel  ay  to  remote. 

750  LOCAL  LOCKOUT  5  !  Local  lockout  for  DVM  and  relay. 

7S0  OUTPUT  502  ;MF0R2T0X"  !  DVM  to  DCV,2V  range,  and  cont.  on  talk. 

770  OUTPUT  507  **A12M  !  Connect  all  A1&A2  terminals  to  C1&C2  terminals. 

780  ENTER  502  ;  AS  !  Trial  read  of  DVM. 

790  IJ1=0  !  Set  counter  for  number  of  times  power  turned  on, 

800  IJ2=0  !  Set  counter  for  number  of  times  power  turned  off. 

810  J=1  !  Set  counter  for  reading  probe  B. 

820  NT=NB/NA  !  Set  interval  between  reading  probe  B. 

830  IF  KK>0  THEN  880 

840  SETTIME  0,0  !  Zero  clock  for  easy  interpretation  only. 

850  !  Loop,  read  probe  A,  read  time,  put  values  into  arrays. 

860  NNA= 1  ?  NNB= 1 

870  FOR  1=1  TO  AJOBS  9  ENTER  5C2  ;  AS9  ATIME< I >=TIME  +<DATE  -EXPDATE > -86400  »  A' 
EMP<I)=VAL  < A$ [ 5 1 ) 

880  ATEMP( I )  =PR0BC 1 -ATEMP ( I ) +PR0BC2  !  Evaluate  temp,  from  voltage  reading. 

890  NNA-NNA+1 

900  DISP  "ATEMP  =  "  ,  ATEMP  ( I )  ,  "TIf1E  =  "  ,  ATIME<  I ) 

910  IF  I=J*NT  THEN  930  !  Check  to  see  if  probe  B  is  to  be  read. 

920  GOTO  1010 

930  OUlPUT  507  "Bl"  !  Switch  B1  to  Cl  prior  to  reading  probe  B. 

940  WAIT  200  !  Wait  200ms  for  DVM  to  read  probe  B. 

950  !  Read  probe  B,  read  time,  put  values  into  arrays.  Switch  A1  to  Cl. 

960  ENTER  502  ;  BS*  BT IME ( J> =TIME  +(DATE  -EXPDATE >*86400  9  BTEMP(J>=VAL  <B$[5i> 

*  OUTPUT  507  "A  1 ” 

970  BTEMP ( J )  =PR0BC3-BTEMP < J ) +PR0BC4  !  Evaluate  temp,  from  voltage  reading. 

980  NNB-NNB+1 

990  DISP  "BTEMP=" ,BTEMP< J> , "T IME=" , BT IME( J )  9  DISP  9  DISP  9  DISP  9  KEY  LABEL 
1000  J  =  J+1  !  Increment  counter  on  read  probe  B. 

1010  WAIT  NA  ?  Wait  time  between  reading  proDe  A. 

1020  NtXT  I  !  Increment  loop. 

1030  LOCAL  5 

1040  OUTPUT  702  "B12"  ?  Simply  turn  off  lights  on  relay. 

1050  DISP  "Storing  data  on  disc.  You  will  be  prompted  for  more  inputs  shortly." 

1060  DISP  9  DISP  9  DISP 

1070  WAIT  2000 

1C80  60SUB  DISCPRINT 

1090  GOTO  1400 

M00  CLEAR. 

1110  DISP  "Do  not  press  any  keys  until  the  computer  responds  witn  a  beep  and" 
1120  DISP  displays  'Data  storage  is  complete'.  You  will  then  be  prompted  to" 
1130  DISP  "key  in  new  observation  data." 

1140  WAIT  5000 

I1 50  GOSUB  DISCPRINT 

1160  GOTO  440 

1  1 70  DISCPRINT:  !  Write  all  data  to  disc  for  this  observation  interval. 

1180  PRINT--  1  ;  NNA-1 
1  1  90  FOR  I  =  i  Tu  NNA- 1 
1200  PRINT*  1  ;  ATEMP < I > . AT IME  < I > 


1210  NEXT  I 

1220  PRINT*-  1  ;  NNB- ! 

1230  FOR  1=1  TO  NNB-1 

1240  PRINT*-  1  ;  BTEMP(I)  ,BT1ME<I> 

1250  NEXT  I 

1280  PRINT#  1  ;  IJ1  ?  The  number  of  times  power  turned  on. 

1270  IF  IJ1=0  THEN  1290 

1280  FOR  1=1  TO  IJ1  «  PRINT#  1  ;  PONT<I)  9  NEXT  I 

1290  PRINT#  1  :  IJ2  !  The  number  of  times  power  turned  off. 

1300  IF  IJ2-0  THEN  1320 

1310  FOR  1=1  TO  IJ2  9  PRINT#  1  :  POFFT(I)  9  NEXT  I 
1320  PRINT#  1  ;  TODS  !  Writes  time  of  day  to  disc 
1330  KK=KK+1 
1340  BIN=AJDBS 

1350  NDPTS=NDPTS+2*( NNA+NNB) +1 J 1 +1 J2+ 1 0 

1360  BEEP  6  CLEAR  6  DISP  "Data  storage  is  complete."  9  RETURN 
1370  GOTO  1400 

1380  BEEP  9  IJ1-IJ1+1  9  DISP  "Power  on"  9  PONT  ( IJ1 )  =TIME  9  RETURN 
1390  BEEP  6  IJ2-IJ2+1  8  DISP  "Power  off"  9  POFFTf IJ2)=TIME  €  RETURN 
1400  DISP  "Do  you  wish  to  record  more  data?  ENTER  (Y)ES  or  (N>0." 

1410  INPUT  REPLYS 

1420  IF  REPLYS="Y"  THEN  420 

1430  ASSIGN#  1  TO  * 

1440  DISP  "The  final  disc  management  is  now  being  completed" 

1450  ASSIGN#  1  TO  EXPDATAS6" : D50 1 " 

1460  IA 1  =  1  9  IAT0T=0  9  IB  1  =  1  9  IBTOT  =  0  9  ION  1  =  1  9  IONT=0  9  I0FF1  =  !  9  I0FFT  =  Q 
1470  FOR  J1=1  TO  KK 

1480  READ#  1  ;  ANUM<J1)9  IAT0T=IAT0T+ANUM( J1 ) 

1490  FOR  I-IA1  TO  IATDT  $  READ*  1  ;  ATEMP(I) ,ATIME(I)9  NEXT  I 
1500  READ*  1  ;  BNUM(J1>9  IBTGT= IBTQT+BNUM< J1  ) 

1510  FOR  I  =  IB  1  TO  IBT0T  9  READ*  1  ;  BTEMP(I) ,BTIME(I)9  NEXT  I 

o2U  RlAu#  l  ;  lrtiNiJlis*  lUN  i  =  lDN  l  ^irUN(  J 1  )  9  IF  IP0N(J1>  =  G  THEN  1540 


1530 
1540 
1550 
1560 
1570 
1580 
1590 
1600 
1610 
1620 
1630 
1640 
1650 
1  660 
1670 
1680 
1690 
T  9 


FOR  I =  I  ON  1  TO  ION  I  w  klhu#  i  :  r  G»i  <  f  i  >  t?  is 
READ#  1  ;  IP0FF( J1 )9  I0FFT=I0FFT+IP0FF< J1 ) 
FOR  I  =  I OFF  1  TO  IOFFT  9  READ*  1 
IAl =IA 1 +IATQT  9  IB1  =  IB  1 +IBT0T  9 
NEXT  J1 

READ#  1  ;  TODS 
ASSIGN*  1  TO  * 

PURGE  EXPDATA£d":D501" 

PACK  "  *.  D501" 

NDPTS=8*NDPTS 

CREATE  EXPDATAS&" : D50 1" , 1 .NDPTS 
ASSIGN#  1  TO  EXPDATAS&" : D50 1 " 


i_  A  i  ± 


9  IF  IP0FF ( J1 >=0  THEN  1560 
P0FFT  <1)9  NEXT  I 

ION  1  =  IQN1 +I0NT  9  I0FF 1  =  IOFF 1 +I0FFT 


PRINT*  1  ;  IATOT  9  DISP  "TOTAL  PROBE  A  READINGS  -  " . IAT0T 

FOR  1=1  TO  IATOT  9  PRINT*  1  ;  ATEMP(I) ,ATIME(I>  9  NEXT  I 

PRINT*  1  ;  IBTOT  9  DISP  "TOTAL  PROBE  B  READINGS  «  ",IBT0T 

FOR  1=1  TO  IBTOT  9  PRINT#  1  ;  BTEMP< I > . BT IHE < I ) • 9  NEXT  I 

PRINT*  1  ;  10NT  9  DISP  "MW  POWER  TURNED  ON  ",I0NT,"  TIMES" 

PRINT#  1  ;  PONT ( I >  9  DISP  P0NT<I>  9  NEXT  I 

1  ;  IOFFT  9  DISP  "MW  POWER  TURNED  OFF  ", IOFFT,"  TIMES"  9  FOR  I- 
1  ;  POFFT(I)  9  DISP  POFFT(I)  9  NEXT  1 
Time  of  day", TODS 


1700  PRINT* 

IOFFT  9  PRINT* 

1710  PRINT*  1  ;  TODS  9  DISP 
1720  ASSIGN*  1  TO  - 
1730  DISP  "Sample  number  is 
1740  DISP  "Date  of  exposure 
1750  DISP  "Finished  the  disc 
*760  END 


9  FOR  1  =  1  TO  IOt 


TO 


\  SAMPLES 
is  " .EXPDATES 
management.  GOODBYE 


HAVE  A  NICE  DAY  1 


PROGRAM  FOR  DATA  ANALYSIS,  CURVE  SMOOTHING  AND  GRAPHICS  _  version  A 
]Q  t  ****«#**« * * * **•* *  * ************************* * ***** 

20  !  *  * 

30  »  *  PROGRAM  TO  READ  TEMPERATURE  READINGS  TAKEN  FROM  * 

40  !  *  TWO  NARDA  PROBES  AND  STORED  ON  DISC,  PLOT  DATA  * 

50  »  *  IN  VARIOUS  FORMS  AND  DO  DOSIMETRY  CURVE  FITTING  * 


MO  ! 
50  » 
60  » 
70  ! 
80  ! 


-  WRITTEN  BY  CHRISTOPHER  C  DAVIS,  DECEMBER,  1986 


90  }  *«*«**»**»***,»«»*»««««»*«**********«*******«**«***** 

100  OPTION  BASE  1 

110  MASS  STORAGE  IS  ":D500" 

120  DIM  ATEMP( 1000) ,ATIME< 1000) ,BTEMP< 1000) ,BTIME( 1000) ,PONT< 50) ,POFFT(50) 
130  DIM  ANUM(20 ) , BNUM(20 ) ,IPON(20) ,IPOFF(20) , DUMMYSt  50 ] 

1 40  !  ***  FUNCTION  TO  GIVE  INVERSE  VIDEO  DISPLAYS 
150  DEE  FNIVID$(DUMKY$> 

160  LENGTH=LEN  (DUMMYS) 

170  FOR  11=1  TO  LENGTH 

180  DUMMYSt II , II]-CHRS  <NUM  (DUMMYSt II , II ] >+128) 

190  FNIVID$=DUMMY$ 

200  NEXT  II 
210  FN  END 
220  CLEAR 

230  DISP  FNIVID$( "DATA  DATE'?") 

240  DISP  FNIVIDSCDATE  FORMAT") 

250  DISP  FNIVIDSC "DAY  MONTH, YEAR") 

260  INPUT  EXPDATES 

270  DISP  FNIVIDS( "ENTER  SAMPLE  ID") 

280  INPUT  SAMPLES 
290  EXPDATA$= EXPDA TE5&SAMPLES 
300  ON  ERROR  GOTO  320 
310  GOTO  370 

320  IF  ERRN  =67  THEN  DISP  "YOU  HAVE  ENTERED  THE  WRONG  FILE  NAME" 

330  IF  ERRN  #67  THEN  DISP  "AN  ERROR  HAS  OCCURRED  -  START  AGAIN" 

340  IF  ERRN  #67  THEN  GOTO  10 

350  CAT  ":D501" 

360  DISP  "THESE  ARE  THE  DATA  FILES  -  CHOOSE  THE  ONE  YOU  WANT"  9  GOTO  230 
370  ASSIGN#  1  TO  EXPDATAS&" : D50 1 " 

380  NUMSMO=0  »  NTEST=0  £  SAR=0  9  T0D$="0" 

390  IAT0T  =  0  &  IBTOT  =  0  9  I0NT  =  0  %  I0FFT=0  «  I AT  =  T  6  IB1=1  €  I0N1-1  9  IOFF 1  =  T 

400  READ#  1  :  IATOT  1  READ  NO  OF  PROBE  A  READINGS 


410  FOR  I=IA1  TO  IATOT  £  READ#  1  ;  ATEMP(I) ,ATIME(I)%  NEXT  I  !  PROBE  A  DATA 

420  READ#  I  ;  IBTOT  !  READ  NO  OF  PROBE  B  READINGS 

430  FOR  I=IB1  TO  IBTOT  £  READ#  1  ;  BTEMP(I) ,BTIME(I>£  NEXT  I  !  PROBE  B  DATA 

440  READ#  1  ;  IGNT  !  NO  OF  TIMES  MW  PWR  TURNED  ON 

450  FOR  I=ION1  TO  IONT  9  READ#  1  ;  PONT<I)£  NEXT  I  !  ON  TIMES 

460  READ#  1  ;  IOFFT  !  NO  OF  TIMES  MW  PWR  TURNED  OFF 

470  FOR  I- IOFF 1  TO  IOFFT  9  READ#  1  ;  POFFT<I>&  NEXT  I  !  OFF  TIMES 

480  IAWA1 +IATOT  9  IB  1 =IB1 +IBTOT  9  ION1 = ION  1 +I0NI  9  IOFF 1 = IOFF 1 + IOFFT 

490  ON  ERROR  GOTO  499 

491  READ#  1  ;  TODS 

499  OFF  ERROR  A  ASSIGN#  1  TO  * 

500  !  PURGE  EXPDATAS4" :D501" 

510  DISP  "TOTAL  PROBE  A  READINGS  -  ".IATOT 

520  ?Jt2£  1-1  T0  IAT0T  •  DISP  ATEMP<I),ATIME<I),I  &  NEXT  I 

530  DISP  "TOTAL  PROBE  B  READINGS  -  ".IBTOT 

ten  LoSR..Im1  T0  IBT0T  •  DISP  BTEMP<  I),BTIME(I),I  9  NEXT  I 
550  DISP  "MW  POWER  TURNED  ON  ".IONT,"  TIMES" 

560  DISP  "MW  POWER  TURNED  OFF  ", IOFFT,"  TIMES" 

toS  rnn^T  *'THE  MW  PWR  WAS  TURNED  ON  AT  THE  FOLLOWING  TIMES" 

|f0  FOR  1-1  TO  IONT  9  DISP  PONT(I)  £  NEXT  I 
590  DISP  "Savoie  ID  is  " ,FNIVIDS(SAMPLES > 


600  DISP  "Date  of  exposure  was  " . FNIVIDS( EXPDATES) 

601  IF  TOD$="0"  THEN  GOTO  610 

602  DISP  "Time  of  day  was  " , FNIVIDSt T0D$> 

610  !  **«**«**>•>«***«**>•»«**«*****«*««*«»****«««**«**<*** 

620  !  *  SECTION  FOR  SMOOTHING 

630  !  *  AND  PLOTTING  OF  EXPERIMENTAL  DATA 

640  !  *  * 

650  !  *  WRITTEN  BY  Christopher  C.  Davis,  June  1934.  * 

660  !  *  MODIFIED  DECEMBER  1986 

670  t  * 

680  !  *«**********»**•«***«*«*»***«*******<•«**«**»***<***« 

690  ! 

700  ! 

710  DIM  AA(5) ,SN(21 ,21 > ,R(20) ,HEAT<5) , COOL (5)  !  ARRAYS  USED  IN  CURVE  SM00THIN 
G 

720  NYSLS=”N" 

730  IF  ATIME(IAT0T)>BTIME(13T0T)  THEN  XMAX  =  AT IME ( I A  TOT )  ELSE  XMAX=BTIME( IBTuI 
)  !  SET  X  AXIS  MAX  TO  MAX  OBSERVATION  TIME 
740  XMIN=0  !  INITIALLY  SET  X  ORIGIN  TO  0 
750  JMIN= 1  §  KMIN= 1 

760  DISP  "TEMPERATURE  MEASUREMENT  INTERVAL  WAS",XMAX 
770  YMIN= 1 0' 20  *  YMAX  =  - ( 1 0~20  > 

780  !  «**  FIND  MAXIMA  AND  MINIMA  OF  DATA 
790  FOR  1=1  TO  IATOT 

800  IF  ATEMP<I)>YMAX  THEN  YMAX=ATEMP< I ) 

810  IF  ATEMPdXYMIN  THEN  YMIN  =  ATEMP<  I ) 

820  NEXT  I 

830  FOR  1=1  TO  IBTOT 

840  IF  BTEMP( I)>YMAX  THEN  YMAX=BTEMP( I) 

850  IF  BTEMPUXYMIN  THEN  YMIN  =  BTEMP(  I ) 

860  NEXT  I 

870  NUMPLT=0  $  IF  NTEST*0  THEN  CLEAR 

380  DISP  "DO  YOU  WISH  TO  PLOT  PROBE  A,  PROBE  B,  OR  BOTH?  "&FNIVIDS< "ENTER  1,2 
OR  3") 

890  DISP  "YOU  CAN  ONLY  CALCULATE  SAR  USING  PROBE  A" 

900  INPUT  NPROBE 

910  IF  NPRQBE= 1  THEN  GOTO  930 

920  GOTO  940 

930  FOR  1-1  TO  IATOT  $  BTIME(I)=ATIME(I)  «  BTEMPt  I )  =ATEMP  ( I )  Li  NEXT  I 
940  IF  NTEST=0  THEN  GOTO  960 

950  NYSLS=" Y" 

960  DISP  "ENTER  TIME  REGION  YOU  WISH  TO  ANALYZE" 

970  INPUT  TMIN.TMAX 

980  FOR  1=1  TO  IATOT  !  FIND  RANGE  OF  TEMP  ARRAY  TO  BE  PLOTTED 
990  IF  TMIN>ATIME ( I )  THEN  JMIN-I 
1000  NEXT  I 

1010  FOR  1=1  TO  IBTOT  9  IF  TMIN>BTIME< I )  THEN  KMIN-I 
1020  NEXT  I 

1030  IF  ATIME(JMINXBTIME(KMIN>  THEN  TMIN=ATIME< JMIN)  ELSE  TMIN  =  BTIME(KMIN) 
1040  FOR  1=1  TO  IATOT  %  IF  TMAX>AT IME( I )  THEN  JMAX-I 
1050  NEXT  I 

1060  FOR  1=1  TO  IBTOT  9  IF  TMAX>BTIME< I )  THEN  KMAX-I 
1070  NEXT  I 

1080  IF  ATIMEt JMAX)>BTIME(KMAX>  THEN  TMAX=ATIME< JMAX)  ELSE  TMAX=BT IME(KMAX) 

1090  IF  NTEST-0  THEN  GOTO  1400 

1100  IF  NYSL$*"N"  OR  NYSL$-"NO"  THEN  GOTO  1400 

11  IP  ?  ***  SECTION  FOR  SELECTION  OF  CURVE  REGION  FOR  SMOOTHING 

1120  TEMAX*-1.E20  9  TEMIN-1.E20 

1130  FOR  I -JMIN  TO  JMAX  9  IF  ATEMPt  I XTEMIN  THEN  TEMIN=ATEMP< I ) 

1140  iF  ATEI1P<  1 1-TEMIN  THEN  JTURN-I  !  Find  inflection  point  on  temp  curve 
1151)  NEXT  I 


1160  DISP  "NUMBER  OF  POINTS  ON  COOLING  SLOPE-  ",  JTURN- JMIN  + 1 

1170  DISP  "NUMBER  OF  POINTS  ON  HEATING  SLOPE-  " , JMAX-JTURN+1 

1180  DISP  "IS  THIS  SATISFACTORY?  Y(ES)  OR  N(0>"  9  INPUT  NYSAT$ 

1190  IF  NYSAT$="N"  OR  NYSATS="NO"  THEN  GOTO  960 

1200  IF  NYSL5="N"  OR  NYSLS="NO"  THEN  GOTO  MOO 

1210  CLEAR  9  DISP  "COOLING  REGION  POINTS" ,JMIN , JTURN 

1220  DISP  "HEATING  REGION  POINTS", JTURN, JMAX 

1230  ON  KEY#  1 ."COOLING"  GOTO  1280 

1240  ON  KEY#  2, "HEATING"  GOTO  1340 

1250  DISP  FNIVID$<"SELECT  OPTION")  9  DISP  9  DISP 

1260  KEY  LABEL 

1270  GOTO  1270 

1280  NSM-1  9  DISP  "ENTER  POINTS  RANGE  FOR  COOLING  SLOPE" 

1290  DISP  "ENTER  0,0  FOR  HELP  IN  FINDING  TURNING  POINT  IN  HEATING  CURVE"  9  INPUT 
JMIN, JTURN 

1300  IF  JMIN-0  AND  JTURN-0  THEN  GOSUB  Find 
1310  DISP  "JTURN-", JTURN 

1320  DISP  "ENTER  JMIN. JTURN"  $  INPUT  JMIN, JTURN 
1330  GOTO  1390 

1340  NSM=2  9  DISP  "ENTER  POINTS  RANGE  FOR  HEATING  SLOPE" 

1350  DISP  "ENTER  0,0  FOR  HELP  IN  FINDING  TURNING  POINT  IN  HEATING  CURVE"  6  INPUT 
JTURN, JMAX 

1360  IF  JTURN-0  AND  JMAX-0  THEN  GOSUB  Find 
1370  DISP  "JTURN-", JTURN 

1380  DISP  "ENTER  JTURN, JMAX"  9  INPUT  JTURN, JMAX 

1390  NUMPlT=1  9  GOSUB  Smooth 

1400  GOSUB  Aplot 

1410  ON  KEY#  1, "AGAIN"  GOTO  1530 

1420  ON  KEY#  2, "DOSE"  GOTO  1580 

1430  ON  KEY#  3, "PLOT"  GOTO  1590 

1440  ON  KEY#  4,"SCRDMP"  GOTO  3910 

1450  ON  KEY#  5, "EXIT"  GOTO  3920 

1460  CLEAR  9  KEY  LABEL 

1470  DISP  "SELECT  OPTION  -  "<iFNIVID$("DOSE")A"  ALLOWS  CALCULATION  OF  SAR" 

1480  DISP  "  -  "&FNIVID$("AGAIN" )&"  RE-READS  DATA  FROM  DISC" 

1490  DISP  "  -  "aFNlVID$( "PLOT" )&"  RE-PLOTS  LAST  PLOT  -  USE  WITH  C 

ARE” 

1500  DISP  "  -  "&FNIVID$< "SCRDMP"><£"  DUMPS  GRAPHICS  TO  PRINTER" 

1510  DISP  "  -  "&FNlVIDS("EXIT")i"  TERMINATES  EXECUTION" 

1520  GOTO  1520 

1530  CLEAR 

1540  DISP  "IF  YOU  CONTINUE  DATA  WILL  BE  RE-READ  FROM  DISC" 

1550  DISP  "PRESS  "&FNIVIDS("AGAIN")&"KEY  IF  YOU  WISH  TO  DO  THIS" 

1560  ON  KEY#  1, "AGAIN"  GOTO  370  9  KEY  LABEL 
1570  GOTO  1570 
1580  GOTO  930 

1590  IF  NUMSMO-O  THEN  GOTO  1620 

1600  GOSUB  Aplot  9  GOTO  1520 

1610  GOTO  1410 

1620  DISP  "OPTION  NOT  AVAILABLE"  9  GOTO  370 
1630  Smooth: 

1640  ! . - . - . . 

1650  !  * 

1660  !  -  CURVE  SMOOTHING  SECTION 

1670  !  «  * 

1680  1  . . 

1690  DISP  "CURVE  FITTING  IN  PROGRESS  "AFNIVID* ("PLEASE  WAIT") 

1700  NOMSMO-1 

1710  IF  NSM-2  THEN  GOTO  1750 

1720  N- JTURN- JMIN+ 1 

1730  FOR  I-i  TO  N  9  BTIME< I >-BTIME( JMIN+I- 1 >  9  BTEMP( I ) -BTEMPt  JMIN+ 1- 1 >  9  NEX 


T  I 
1  740 
1750 
1760 
RN+I-1) 
1770 
1780 
1790 
1800 
1810 
1820 
1830 
1840 
1850 
1860 
1870 
1880 
1890 
1900 
1910 
1920 
1930 
1940 
1950 
1960 
1970 
1980 
1990 
Ni ) 


GOTO  1780 
N= JMAX- JTURN+1 

FOR  1  =  1  TO  JMAX- JTURN+ 1  9  BTIME< I ) =BT IME( JTURN+I- 1 )  9  BTEMP< I )  =BTEtiP(  JTU 
NEXT  I 

FOR  1=1  TO  5 

R(I)=°  $  AA( I ) =0  9  FOR  J=1  TO  21  9  SN(I,J)=0  9  NEXT  J  9  NEXT  I 
DISP  WHAT  ORDER  OF  POLYNOMIAL  DO  YOU  WISH  TO  FIT  -  1 'LINEAR ... " 

INPUT  NS 

FOR  J=1  TO  NS+1 

FOR  1=1  TO  N 

IF  J=NS+ 1  AND  BTIME( I)=0  THEN  GOTO  1860 

R( J)  =  BTEMP (I)*BTIME(I)'(NS+1-J)+R(J)  9  GOTO  1870 

R< J)=BTEMP( I)+R( J) 

NEXT  I 
NEXT  J 

FOR  K=1  TO  NS+1  9  FOR  J=1  TO  NS+1  9  FOR  1=1  TO  N 
IF  2*NS- J-K+2=0  AND  BTIME<I)=0  THEN  GOTO  1920 
SN( J,K) =BTIME< I ) ' <2*NS- J-K+2) +SN< J,K)  9  GOTO  1930 
SN  <  J  t  K ) *  1 +SN  t  J ,  K  > 

NEXT  I 
NEXT  J 
NEXT  K 
GDSUB  Invt 

FOR  1  =  1  TO  NS+1  at  FOR  J=l  TO  NS+1 

AA(NS+2-I)=SN(NS+2-I.J)*R(J)+AA(NS+2-I)  9  NEXT  J  9  NEXT  I 

IF  NSM= 1  THEN  DX=ATIME< JTURN) -ATIME( JMIN)  ELSE  DX=ATIHE< JMAX)-ATIME( JTUR 


2000  IF  NSf1=  1  THEN  XSTART  =  ATIME( JMIN)  ELSE  XSTART  =AT IME< JTURN) 

2010  FOR  K=1  TO  200  «  I=K+800  BTIME < I ) = <K- 1 )*DX/200+XSTART  «  BTEMP(I)=0  e  F 
OR  J=1  TO  NS+1 

21)20  IF  NS+ 1  - J  =  0  AND  BTIME<I)  =  0  THEN  GOTO  2050 

2030  BTEMP(I)=AA(J)»BTIME<I)' ( NS+ 1  - J) +BTEMP( I ) 

2040  GOTO  2060 

2050  BTEMP(I)=AA( J)+BTEMP(I) 

2060  NEXT  J 

2070  NEXT  K 

2080  RETURN 

2090  Invt : 

2100  ! 

2110  »  MATRIX  INVERSION  BY  GAUSS-JORDAN  ELIMINATION  WITH  FULL  PIVOTING 

2120  !  SQUARE  MATRIX  OF  ORDER  NP 

2130  ! 

2140  DIM  KP(5) , JP(5) 

2150  NP=NS+ 1 

2160  NCOL-NP 

2170  ! 


2180 
2190 
2200 
2210 
2220 
2230 
2240 
2250 
2260 
227  0 
2280 
2290 
2300 
2310 


!  INITIALIZE  PERMUTATION  VECTORS 

* 

FOR  1=1  TO  NP  9  KP( I ) *  I  9  JP(I)=I  «  NEXT  I 

» 

!  INVERT  MATRIX 

» 

FOR  IR- 1  TO  NP 

f 

*  FIND  THE  PIVOT  ELEMENT 

t 

PIVOT  =  0 

FOR  1  -  IR  TO  NP  9  FOR  J=IR  TO  NP  9  K-KPU)  9  L  =  JP(J) 
TEST  =ABS  <SN<K.L))-ABS  (PIVOT) 

IF  TEST <  =  0  THEN  GOTO  2370 


2320  IF  TEST>0  T!'FN  GOTO  2340 

2330  GOTO  2370 

2340  IK«I 

2350  JK-J 

2360  PIVOT=SN(K ,L) 

2370  NEXT  J  $  NEXT  I 

2380  ! 

2390  !  THE  PIVOT  ELEMENT  IS  SN(K,L> 

2400  ! 

2410  !  UPDATE  PIVOT  VECTORS 

2420  ! 

2430  K=KP<  IK)  S  KP(IK>=KP<IR>  *  KP<IR>=<  C«  L=JP<JK)  £  JP<  JK)=JP(  IR)  JP(IR) 
L 

2440  ! 

2450  !  COMPUTE  NEW  ELEMENTS  OF  PIVOTAL  ROW 

2460  ! 

2470  FOR  J=1  TO  NCOL  d  SN< K , J ) =SN< K , J) /PIVOT  $  NEXT  J 
2480  * 

2490  1  COMPUTE  REMAINING  ELEMENTS  OF  IR' TH  STEP 

2500  ! 

2510  FOR  1=1  TO  NP 

2520  IF  I-K=0  THEN  GOTO  2590 

2530  AIL=SN( I , L ) 

2540  SN(I,L)=-(AIL/PIVOT) 

2550  FOR  J=1  TO  NCOL 

2560  IF  J-L=0  THEN  GOTO  2580 

2570  SN< I , J>=SN< I , J)-AIL*SM(K , J) 

2580  NEXT  J 

2590  NEXT  I 

2600  1 
2610  ! 

2620  !  INVERT  PIVOT  ELEMENT 

2630  ! 

2640  SN(K . L ) *  1 /PIVOT  *  NEXT  IR 

2650  RETURN 

2660  Aplot : 

2670  !  . . . 

2680  !  * 

2690  !  *  PLOTTING  SUBROUTINE  * 

2700  i  . . . 

2710  ! 

2720  NLT=1 

2730  CLEAR  &  IF  NUMPLT=0  THEN  GOTO  2930 
2740  ALPHA  $  DISP  "PARAMETERS  OF  FIT"  $  DISP 
2750  FOR  1  =  1  TO  NP  DISP  AA(I)  $  NEXT  I 
2760  IF  NSM=2  THEN  GOTO  2790 
2770  FOR  1=1  TO  NP  $  COOL ( I ) = AA( I )  &  NEXT  I 
2780  GOTO  2800 

2790  FOR  1=1  TO  NP  @  HEAT < I ) =AA( I >  NEXT  I 
2800  CHI2-0 

2810  IF  NTEST<2  THEN  GOTO  2850 
2820  SAR-4 .2* (HEAT ( 1  )-COOL<  1  >>**1000 

2830  IF  NS*2  THEN  SAR=2*SAR*AT IME ( JTURN) +4 . 2* (HEAT (2) -C00L<2 ) ) 

2840  DISP  FNIVIDSC'SAR  IS" ) . SAR , "W/k g" 

2850  FOR  1=1  TO  N  ^  Y2=0  $  FOR  J-1  TO  NS*1 
2860  IF  NS+1-J-0  AND  BT1ME<I)=0  THEN  GOTO  2860 
2b 70  Y2'  =  AA(JJ-BTIME(I)  (NS+1-J)  +  Y2  9  GOTO  2890 
2880  Y2*AA( J)+Y2 
2890  NEXT  J 

2900  CHI2=(BTEMP(I>-Y2) '2/Y2+CHI2 
2910  NEXT  I 


DISP  "OBSERVED  CHI2=  ",CHI2,"  WITH  ".N+NS+1,"  D.F." 

DISP  "DO  YOU  WANT  A  HARD  COPY  PLOT?  Y(ES)  OR  N(O)” 

INPUT  NYSS>  IF  NUMSMO=  1  THEN  GOTO  3040 
DISP  "TIME  WILL  BE  PLOTTED  FROM  " , XMIN , "TO" , XMAX 
DISP  "TEMPERATURE  WILL  BE  PLOTTED  FROM" , YMIN,"TO’  » YMAX 
DISP  "IF  YOU  ARE  SATISFIED  PRESS  RETURN" 

INPUT  NYLIMS 

IF  NYLIM$=""  THEN  GOTO  3040 

DISP  "ENTER  TIME  RANGE  TO  BE  PLOTTED,  TMIN.TMAX" 

INPUT  XMIN, XMAX 

DISP  "ENTER  TEMPERATURE  RANGE  TO  BE  PLOTTED" 

INPUT  YMIN.YMAX 

IF  NTEST>=  1  THEN  GOTO  3340 

CLEAR 

DISP  "CALCULATING  PLOT  SCALING  PARAMETERS  "&FNIVIDS( "PLEASE  WAIT") 
AXMAX=INT  (LGT  (XMAX)) 

AYMAX=INT  (LGT  (YMAX)) 

XF=1 0' AXMAX  6  YF=10  AYMAX  «  DELX=.01  e  DELY=.01 


2930 

2940 

2950 

2960 

2970 

2980 

2990 

3000 

3010 

3020 

3030 

3040 

3050 

3060 

3070 

3080 

3090 

3100 

3110 

3120 

3130 

3140 

3150 

3160 

3170 

3180 

3190 

3200 

3210 

3220 

3230 

3240 

3250 

3260 

3270 

3280 

3290 

3300 

3310 

3320 

3330 

3340 

3350 

3360 

3370 

3380 

3390 

3400 

3410 

3420 

3430 

3440 

3450 

3460 

3470 

3480 

3490 

3500 

3510 

3520 


1  D.F." 
N(0)" 


XPWRS=VAL$  (AXMAX)  1 
!  XMIN=0  G  YMIN=0 
IF  (XMAX- XMIN >/XF>=  .1  THEN  DELX=.02 
IF  (XMAX-XMIN) /XF>=  .3  THEN  DELX=.05 
IF  ( XMAX-XMIN) /XF>=  .5  THEN  DELX=.1 
IF  ( XMAX-XMIN) /XF>=  1  THEN  DELX=.2 
IF  (XMAX -XMIN )/XF>=  3  THEN  DELX=.5 
IF  (XMAX-XMIN)/XF>=  5  THEN  DELX=1 
IF  ( YMAX- YMIN) / YF>=  .1  THEN  DELY=.02 
IF  (YMAX-YMIN)/YF>=  .3  THEN  DELY-. 05 
IF  ( YMAX- YMIN)/YF>=  .5  THEN  DELY=  .  1 
IF  ( YMAX-YMIN)/YF>=  1  THEN  DF.LY-.2 
IF  (YMAX-YMIN)/YF>=  3  THEN  DELY= . 5 
IF  ( YMAX-YMIN)/YF>=  5  THEN  DELY= 1 
DISP  "THE  X-AXIS  TICK  SPACING  WILL  BE" , DELX 
DISP  "IF  YOU  ARE  SATISFIED  PRESS"&FNIVIDS("ENDLINE" ) 

INPUT  NTKS'S  IF  NTKS  =  ""  THEN  GOTO  3280 

DISP  "ENTER  X-AXIS  TICK  SPACING"  C  INPUT  DELX 

DISP  "THE  Y-AXIS  TICK  SPACING  WILL  BE" , DEL Y 

DISP  "IF  YOU  ARE  SATISFIED  PRESS"AFN1VID£("ENDLINE" ) 

INPUT  NTK$@  IF  NTK$=""  THEN  GOTO  3320 

DISP  "ENTER  Y-AXIS  TICK  SPACING"  *  INPUT  DELY 

XNUM=CEIL  (( XMAX-XMIN) /(XF*DELX)) 

YNUM-CEIL  ( ( YMAX-YMIN) /( YF*DELY) ) 

IF  NY$="N"  THEN  GOTO  3380 
PLOTTER  IS  505 

GRAPH  %  IF  NUMPLT-0  THEN  GCLEAR  $  DEG 
GOTO  3420 
PLOTTER  IS  1 
GRAPH 

♦  IF  NUHPLT#0  AND  NTEST>0  THEN  GOTO  3110 

IF  NUMPLT-0  THEN  GCLEAR 

LOCATE  30, RATIO  *100-10,20,95 

!  IF  NUMPLT#0  AND  NTEST>0  THEN  GOTO  3110 

SCALE  XMIN/XF , XMIN/XF+XNUH*DELX , YMIN/YF , YMIN/ YF+YNUM* DELY 

IF  NTEST>0  AND  NUMPLT#0  THEN  GOTO  3670 

FXD  2 

LGRID  DELX, DELY, XMIN/XF, YMIN/YF 
AXS=(XMAX-XMIN)/(XF*10) 

AYSM  YMAX- YMIN) /(YF-10) 

MOVE  (XMAX*XMIN)/(2*XF) , YMIN/YF -2*AY°  «  LORG  b 
LABEL  USING  "K"  ;  "TIME  (SECS)  x  lE’LiXPWRS 
MOVE  XMIN/XF-2*AXS, ( YMAX+YMIN) /(2*YT  ) 


<»  YPWR$  =  VAL5  (AYMAX) 


3S30  LORG  4  3  DEG  @  LDIR  90 

3540  LABEL  USING  ,,KM  ;  "TEMP .  (DEGS  C)  x  1E"&YPWR$ 

3550  LDIR  0 

3560  IF  NPR0BE=2  THEN  GOTO  3630 
3570  LINE  TYPE  1 

3580  MOVE  ATIME< JMIN)/XF , ATEMP< JMIN>/YF 
3590  FOR  I-JMIN  TO  JMAX 
3600  DRAW  AT IME< I >/XF , ATEMP( I>/YF 
3610  NEXT  I 

3620  IF  NPROBE= 1  THEN  GOTO  3660 
3630  LINE  TYPE  5 

3640  MOVE  BTIME(KMIN>/XF .BTEMP <KMIN>/YF 

3650  FOR  I=KMIN  TO  KMAX  9  DRAW  BT IME< I ) /XF , BTEMP < I ) /YF  9  NEXT  I 
3660  IF  NYSL$="N"  OR  NYSL$=”NO"  THEN  GOTO  3720 
3670  LINE  TYPE  4 

3680  MOVE  BTI ME (801 )/XF, BTEMP <801 )/YF 

3690  FOR  1=801  TO  1000 

3700  DRAW  BTI ME (I>/XF, BTEMP <I)/YF 

3710  NEXT  I 

3720  PEN  UP 

3730  MOVE  XMIN/XF , YMIN/YF 

3740  IF  SAR=0  THEN  GOTO  3780 

3750  MOVE  XMIN/XF+4 .5*AXS, YMAX/YF+ . 1«AYS 

3760  LORG  4 

3770  LABEL  ”SAR",SAR,”W/kg" 

3780  NTEST=NTEST+1 

3785  DISP  "PRESS  "<iFNIVID$("CONTINUE")A"TO  CALCULATE  SAR" 

3790  PAUSE 
3800  KEY  LABEL 
3810  RETURN 

3820  Find:  !  *—  SUBROUTINE  TO  FIND  POINT  WHERE  HEATING  STARTS 
3830  DISP  "ENTER  A  TEST  RANGE  J1 , J2"  9  INPUT  J1,J2 
3840  DISP  "TIME", "TEMP".  "DT/Dt" 

3850  FOR  I=J1  TO  J2-1 

3860  DISP  ATIME(I) ,ATEMP(I) , < ATEMP< 1+ 1 ) -ATEMPC I ) ) /< ATIME( 1+ t > -ATIME < I ) > 
3870  NEXT  I 

3380  DISP  "TRY  ANOTHER  RANGE  -  Y(ES)  OR  N(0>"  9  INPUT  NYR$ 

3890  IF  NYR$="Y"  OR  NYRS=" YES"  THEN  GOTO  383C 
3900  DISP  "ENTER  JTURN"  9  INPUT  JTURNp  RETURN 
3910  CHAIN  "DUMP" 

3920  END 


PROGRAM  FOR  DATA  ANALYSIS,  CURVE  SMOOTHING  AND  GRAPHICS  -  version  B 


30  !  *  PROGRAM  TO  READ  TEMPERATURE  READINGS  TAKEN  FROM  * 
40  !  *  TWO  NARDA  PROBES  AND  STORED  ON  DISC.  PLOT  DATA  * 
50  !  *  IN  VARIOUS  FORMS  AND  DO  DOSIMETRY  CURVE  FITTING  * 
60  !  * 


70  !  *  WRITTEN  BY  CHRISTOPHER  C  DAVIS.  DECEMBER,  1986  * 

80  !  *  * 

90  !  *•**■****"**••#**«■••**«** * »*»«•»«*»■««*•**«■* «« * ***«•»•••*•« 

100  OPTION  BASE  1 

110  MASS  STORAGE  IS  ":D500" 

120  DIM  ATEMP <  1 000 ). ATIME( 1 000 ) . BTEMP ( 1 000 ) . BTIMEU  000 ) , PONT < 50 ) , POFFT < 50 ) 
130  DIM  ANUM( 20 )  , BNUM( 20  >  ,IPON(20)  ,IPOFF(20)  , DUMMYS! 50  ] 

MO  !  ***  FUNCTION  TO  GIVE  INVERSE  VIDEO  DISPLAYS 
150  DEF  FNIVIDS(DUMMYS) 

160  LENGTH-LEN  (DUMMYS) 

170  FOR  11=1  TO  LENGTH 

180  DUMMYS [ 1 1 , II ] =  CHRS  (NUM  ( DUMMYS [ II , II ] >  +  128) 

190  FN I VIDS= DUMMYS 
200  NEXT  II 
210  FN  END 
220  CLEAR 


230 
240 
250 
260 
270 
280 
290 
3u0 
3  1  0 
320 
330 

34U 

350 

360 

370 

3B0 

390 

4C0 

410 

420 


DISP  FNIVIDSC'DATA  DATE'?") 

D1SP  FNIVIDS( "DATE  FORMAT") 

DISP  FNIVIDSC "DAY  MONTH, YEAR") 
INPUT  EXPDATES 

DISP  FNIVIDSCENTER  SAMPLE  ID") 
INPUT  SAMPLES 
EXPOATA$=EXPDATES&SAMPLES 
ASSIGN#  1  TO  EXPDATAS&" :D501 " 
NUMSMO-O  *  NTEST  =  0  >«  SAR  =  0 
IATOT  =  0  I B T 0 T  =  0  »  IONT  =  0  * 

IATOT  !  READ  NO  OF 
TO  IATOT  •  READ#  1 
IBTOT  !  READ  NO  OF 
TO  IBTOT  $  READ#  1 
IONT  !  NO  OF  TIMES 
TO  IONT  «  READ#  1 
IOFFT  !  NO  OF  TIMES 


RrAD#  1  ; 

FOR  I = IA 1 
READ#  1  ; 

FOR  I-IB1 
READ#  1  ; 

FOR  1  =  ION  1 
READ#  1 

FOR  I  *  IOFF 1  TO  IOFFT  v*  READ#  1 
IA1 = IA1 + IATOT  ft  IB  I  =  I B 1 +IBT0T  & 
ASSIGN#  1  TO  * 


IOFFT-U  ft  I A 1  =  1  9  IB  1  =  1  ft  ION  1  =  I  ft  IOFF  1  =  I 
PROBE  A  READINGS 

:  ATEMP(I)  .ATIME(I)ft  NEXT  I  !  PROBE  A  DATA 
PROBE  B  READINGS 

;  BTEMP(I)  .  BTIME<  I  >«•  NEXT  I  !  PROBE  B  DATA 
MW  PWR  TURNED  ON 
;  PONT <  I  )ft  NEXT  I  !  ON  TIMES 
MW  PWR  TURNED  OFF 
;  POFFT ( I  >ft  NEXT  I  !  OFF  TIMES 
ION1-ION1+IONT  ft  IOFF  1  =  IOFF  1  +I0FFT 


430  !  PURGE  EXPDATAS&" :D5D1 " 

440  DISP  "TOTAL  PROBE  A  READINGS  =  ".IATOT 

450  !  FOR  1  =  1  TO  IATOT  *  DISP  ATEMP ( I )  ,  ATIME  ( I  > ,  I  ft  NEXT  I 

460  DISP  "TOTAL  PROBE  B  READINGS  -  ", IBTOT 

470  !  FOR  1  =  1  TO  IBTOT  ft  DISP  BTEMP ( I )  ,  BT IME < I >  ,  I  ft  NEXT  I 

480  DISP  "MW  POWER  TURNED  ON  ",IONT,"  TIMES" 

490  DISP  "MW  POWER  TURNED  OFF  ".IOFFT,"  TIMES" 

500  DISP  "THE  MW  PWR  WAS  TURNED  ON  AT  THE  FOLLOWING  TIMES" 
510  FOR  1-1  TO  IONT  ft  DISP  PONT(I)  ft  NEXT  I 
520  DISP  "Sample  ID  is  " ,FNIVID$(SAMPLE5> 

530  DISP  "Date  of  exposure  uas  " ,FNIVID$(EXPDATE$ ) 

540  *  . - . 

-•90  '  *  SECTION  FOR  SMOOTHING 

9bC  !  -  AND  PLOTTING  OF  EXPERIMENTAL  DATA 

5  7  0  1  •  „ 


56U  ' 

590  ! 

600 


«RITTEN  BY  Christopher  C.  Davis. 
MODIFIED  DECEMBER  1966 


June  1984. 


BIO  !  »»»***»*«•*****«»***«*****»*****»******»*«**«**«■•« 

620  * 

630  ! 

640  DIM  AA(5> ,SN<21 ,21 > ,R(20) ,HEAT<5> ,C00L<5>  J  ARRAYS  USED  IN  CURVE  SMOOTHIH 
650  NYSLS="N" 

660  IF  AT]ME<IATOT)>BTIME(IBTOT>  THEN  XMAX=ATIME(IATOT)  ELSE  XMAX=BTIME(IBTQT 
)  !  SET  X  AXIS  MAX  TO  MAX  OBSERVATION  TIME 
670  XMIN-0  !  INITIALLY  SET  X  ORIGIN  TO  0 
680  JMIN-1  6  KMIN= 1 

690  DISP  "TEMPERATURE  MEASUREMENT  INTERVAL  NAS" , XMAX 
700  YMIN-10'20  Q  YMAX=- ( 1 0~20 > 

710  !  FIND  MAXIMA  AND  MINIMA  OF  DATA 
720  FOR  1-1  TO  IATDT 

730  IF  ATEMP( I )>YMAX  THEN  YMAX=ATEMP( I ) 

740  IF  ATEMP(IXYMIN  THEN  YMIN=ATEMP<  I ) 

750  NEXT  I 

760  FOR  1-1  TO  IBTOT 

770  IF  BTEMP( I ) > YMAX  THEN  YMAX=BTEMP< I > 

780  IF  BTEMPdXYMIN  THEN  YMIN  =  BTEMP(I ) 

790  NEXT  I 

800  NUMPLT-0  $  IF  NTEST*0  THEN  CLEAR 

810  DISP  "DO  YOU  NISH  TO  PLOT  PROBE  A,  PROBE  B,  OR  BOTH?  "5FNIVIDS< "ENTER  1,2 
OR  3" ) 

320  DISP  "YOU  CAN  ONLY  CALCULATE  SAR  USING  PROBE  A" 

830  INPUT  NPROBE 

840  IF  NPROBE = 1  THEN  GOTO  860 

850  GOTO  870 

860  FOR  1-1  TO  IATOT  «  BTIMEd ) *HTIME( I )  v  BTEMP(I)=ATEMPCI>  «  NEXT  I 
870  IF  NTEST=0  THEN  GOTO  890 

880  NYSL$="Y" 

390  DISP  "ENTER  TIME  REGION  YOU  WISH  TO  ANALYZE" 

900  INPUT  TMIN.TMAX 

910  FOR  1-1  TO  IATOT  *  FIND  RANGE  OF  TEMP  ARRAY  TO  BE  PLOTTED 
920  IF  TMIN>ATIME ( I )  THEN  JMIN-I 

930  NEXT  I 

940  FOR  1-1  TO  IBTOT  $  IF  TMIN>BT I  ME C I >  THEN  KMIN-I 

950  NEXT  I 

960  IF  ATIME(  JMINXBTIME(KMIN)  THEN  THIN=ATIME( JMIN)  ELSE  TMIN-BTIME(KMIN) 

970  FOR  1-1  TO  IATOT  *  IF  TMAX>AT IME ( 1 )  THEN  JMAX-I 

980  NEXT  I 

990  FOR  1-1  TO  IBTOT  «  IF  TMAX>BTIME<I)  THEN  KMAX=I 

1000  NEXT  I 

1010  IF  ATIME<JMAX>>BTIME<KMAX>  THEN  TMAX=ATIME( JMAX)  ELSE  TMAX-BTIME I.KMAX > 

1020  IF  NTEST=0  THEN  GOTO  1250 

1030  IF  NYSLS="N"  OR  NYSL$-"NO"  THEN  GOTO  1250 

1040  !  ***  SECTION  FOR  SELECTION  OF  CURVE  REGION  FOR  SMOOTHING 

1050  TEMAX- -  1 . E20  »  TEMIN-1.E20 

1060  FOR  I -JMIN  TO  JMAX  $  IF  ATEMP< I XTEMIN  THEN  TEMIN-ATEMP (I ) 

1070  IF  ATEMP( I ) - TEM1N  THEN  JTURN-I  !  Find  inf  lection  point  on  temp  curve 
1080  NEXT  I 

1090  DISP  "NUMBER  OF  POINTS  ON  COOLING  SLOPE-  " , JTURN- JHIN+ 1 

1100  DISP  "NUMBER  OF  POINTS  ON  HEATING  SLOPE-  "  .  JMAX- JTURN-*- 1 

1110  DISP  "IS  THIS  SATISFACTORY?  Y<ES>  OR  N<0>"  «  INPUT  NYSATS 

1120  IF  NYSAT$-"N"  OR  NYSATS»"NO"  THEN  GOTO  890 

1130  IF  NYSL$-"N"  OR  NYSLS-"NO"  THEN  GOTO  1250 

1140  CLEAR  ft  DISP  "COOLING  REGION  POINTS" , JMIN . JTURN 

1150  DISP  "HEATING  REGION  POINTS" , JTURN , JMAX 

1160  ON  KEY*  1, "COOLING"  GOTO  1210 

1170  ON  KEY*  2  .  "HEATING"  GOTO  1230 

118b  DISP  FNIVIDSC’SELECT  OPTION")  DISP  *  DISP 

I  1 90  KEY  LABEL 


1200  GOTO  1200 

(210  N$M= I  $  DISP 

T  JMIN , JTURN 
1220  GOTO  1240 

1230  NSM-2  «  DISP 

T  JTURN, JMAX 


“ENTER  POINTS  RANGE  TOR  COOLING  SLOPE"  9  DISP  &  DISP  S>  INPU 


"ENTER  POINTS  RANGE  FOR  HEATING  SLOPE"  9  DISP  9  DISP  s  INP'J 


1240 

NUMPLT-1 

9  GOSUB  Smooth 

1250 

GOSUB  Aplot 

1260 

ON  KEY# 

1 ."AGAIN"  GOTO  1380 

1270 

ON  KEY# 

2.  "DOSE"  GOTO  1430 

1280 

ON  KEY# 

3, "PLOT"  GOTO  1440 

1290 

ON  KEY# 

4  ,  "SCRDMP"  GOTO  3530 

1300 

ON  KEY# 

5, "EXIT"  GOTO  3540 

1310 

CLEAR  § 

KEY  LABEL 

1320 

DISP  "SELECT  OPTION  -  "4FNIVIDS><  "DOSE"  >d"  i 

1330 

DISP  " 

-  "4FNIV IDS ("AGAIN" >4" 

1340 

DISP  " 

-  "4FNIV IDS< "PLOT" >4" 

ARE" 

1350 

DISP  " 

-  "4FNIVIDS("SCRDMP")4 

1360 

DISP  " 

-  "4FNIVIDS("EXIT")4" 

1370 
1380 
1390 
1400 
1410 
1420 
1430 
!  4^0 
1450 
1460 
1470 
1480 
I4S0 
1500 
1510 
1520 
1530 
1540 
1550 
1  5dU 
1570 
1580 
T  I 
1590 
1600 
1610 
RN+I-1 ) 
1620 
1  630 
1640 
1650 
1660 
1670 
1680 
1690 
1700 
1710 
1  72U 
1  730 
I  74U 
1  750 
1  760 


GOTO  1370 
CLEAR 

DISP  "IF  YOU  CONTINUE  DATA  WILL  BE  RE-READ 
DISP  "PRESS  "4FNIVID£( "AGAIN") 4"KEY  IF  YOU 
ON  KEY#  1, "AGAIN"  GOTO  300  9  KEY  LABEl 
GOTO  1420 
GOTO  860 

IF  NUMSMO=0  THEN  GOTO  1470 


ALLOWS  CALCULATION  OF  SAR" 
RE-READS  DATA  FROM  DISC" 
RE-PLOTS  LAST  PLOT  -  USE  WITH 

"  DUMPS  GRAPHICS  TO  PRINTER" 
TERMINATES  EXECUTION" 


FROM 

WISH 


DISC" 
TO  DO 


THIS" 


GOSUB  Ap lot 
GOTO  1260 

DISP  "OPTION 
Smooth : 

I  **«**«*-***i 


&  GOTO  1370 

NOT  AVAILABLE"  9  GOTO  300 


CURVE  SMOOTHING  SECTION 


•  * 

*  • 

!  * 
f 

DISP  "CURVE  FITTING  IN  PROGRESS  "4FNIVIDS("PLEASE  WAIT") 

NUMSMG- 1 

IF  NSM  =  2  THEN  GOTO  1600 
N= JTURN-JMIN+ 1 

FOR  1-1  TO  N  9  BTIME<  I)-BTIME<  JMIN+I-1  >  *  &TEMP<  I)=BTEMP<  JHIN+I-1  >  e  NEX 


GOTO  1630 

N= JMAX- JTURN+1 

FOR  1-1  TO  JMAX- JTURN+ 1 


9  BTIME<I>*BTIME(JTURN+I-1  )  9  BTEMP  ( I  >  -  BTEI1P  ( JTU 


C 


NEXT  I 
FOR  1-1  TO 
R<I)-0  9  AA< I ) -0 
DISP  "WHAT  ORDER 
INPUT  NS 
FOR  J-1  TO  NS  +  1 
FOR  1-1  TO  N 

IF  J=NS  + 1  AND  BTlME<I)-0  THEN  GOTO  1710 
R< J)-BTEMP( 1 >*BTIHE< I)'  (NS  +  1 -J)+R( J)  <*  GOTO 
R< J)=BTEMP( I)+R( J) 

NEXT  I 
NEXT  J 

FOR  K - 1  TO  NS+1  .4  FOR  J- 1  TO  NS+1  4  FOR  1=1 
IF  2*NS- J-K+2-0  AND  BTIME(I)-U  THEN  GbTu  1/ 
SN<J,K)-BTIME<I>  <  2*NS-  J-k  +3  >  +  SN(  J  ,  K  >  GOT; 


9  FOR  J-l  TO  21  9  SN< I , J)-0  9  NEXT 
OF  POLYNOMIAL  DO  YOU  WISH  TO  FIT  - 


1720 


§  NEXT 
■LINEhk 


TO  N 
'u 

I  1780 


,7  JU 
I  780 
1790 
1  800 
1810 
1820 
1830 
1890 
N) 

1850 
1880 
Oft  J=l 
1870 
1  88  U 
1390 
1900 
1910 
1920 
1930 
1990 
1950 
I960 
1970 
1980 
1990 
2000 
2010 
2020 
2030 
2U9  0 
2050 
2Ub0 
2070 
2U80 
2090 

21  UO 
2110 
2120 
2130 
2190 

2  I  DU 
2160 
2170 
2180 
2190 
2200 
2210 
2220 
2230 
2290 
2250 
2260 
2270 
2280 

2290 

-  Oil  J 

22  I  U 
2-20 
2  o  2  0 

2  2  90 


L 


NEXT  I 
NEXT  J 
NEXT  K 
GOSUB  Invt 

FOR  1=1  TO  NS+1  ft  FOR  J=1  TO  NS  +  1 

AA(NS+2-I>=SN(NS+2-I , J>*R(J>+AA(NS+2-I>  ft  NEXT  J  ft  NEXT  1 

IF  NSM=  1  THEN  DX  =  AT  If1E(  J  TURN  > -A I  THE  ( JtllN)  ELSE  DX=AT  I  ME  C  JMAX ) -AT  DIE  iJTUR 

IF  NSM=  1  THEN  XSTART=ATIME< JMIN)  ELSE  XSTART=ATIME( JTURN) 

FOR  K= 1  TO  200  ft  I=K+800  ft  BT IME ( I ) = ( K- 1 ) *DX/200+XSTART  ft  BT£MP(I)=G  ft  F 

TO  NS+1 

IF  NS  + 1  - J  =  0  AND  BTIME< I) =0  THEN  GOTO  1900 
BTEMP<  I>  =  AA(  J>*6TIME<1) '  (NS+1  -  J)  +tf  i  t.MP(  j  > 

GOTO  1910 

BTEMP ( I >  =  AA<  J) +BTEMP( I ) 

NEXT  J 
NEXT  K 
RETURN 
Invt : 

MATRIX  INVERSION  BY  GAUSS-JORDAN  ELIMINATION  WITH  FULL  PIVOTING 
SQUARE  MATRIX  OF  ORDER  NP 

DIM  KP< 5) , JP(5) 

NP=NS+1 

NCOL=NP 

INITIALIZE  PERMUTATION  VECTORS 

FUR  1*1  TO  NP  £  KP< I ) =1  ft  JP(I)=I  ft  NEXT  I 

INVERT  MATRIX 

FOR  IR= 1  TO  NP 

FIND  THE  PIVOT  ELEMENT 

PIVQT=0 

FDR  I-IR  TO  NP  ft  FOR  J=IR  TO  NP  @  K=KP(I)  ft  L=JP(J) 
lESl=ABS  <SN(K  ,  L  )  7-AB2  (PJ.VUI> 

IF  TEST<=  0  THEN  GOTO  2220 
IF  TEST>0  THEN  GOTO  2190 
GOTO  2220 
IK- 1 
JK=J 

P I VQT  =  SN( K  ,  L  ) 

NEXT  J  ft  NEXT  I 

THE  PIVOT  ELEMENT  IS  SN<K,L> 

UPDATE  PIVOT  VECTORS 

K'KP(IK)  ft  KP(IK)-KP(IR)  ft  KP(IR)-K  ft  L=JP(JK>  ft  JP(jK)=jP(IR)  ft  jP<IR>= 


•  COMPUTE  NEW  ELEMENTS  OF  PIVOTAL  ROW 
(  FGR  >1  TO  NCOL  ft  SN(K,J)-SN(K, J)/PIVGT  ft  NEXT  J 
1  COMPUTE  REMAINING  ELEMENTS  OF  IR'TH  STEP 


L.  o  Cj  \J 

2j  7  0 

2  39  C 
2390 

3  4  u  0 
2410 
2  420 
2430 
24  40 
2450 
2460 
2470 

2480 

2-490 

250G 

2s;o 

2520 

2530 

2540 

2550 

2560 


r  o*  I--4  lb  Nl^ 

if  i-k-0  then  goto  2440 
ail=sn( 1 ,  l  > 

SNi  I ,L  >  =  -  (AlL/PIVGT  > 

F  OR  j- i  10  NCOL 

1-  j-  l  =  0  then  goto  2430 

5N<I  ,  J)-SN(I  ,  J)-AIL*$N(K,  J) 
NEXT  J 
NEXT  I 

t 

f 

i  invert  pivot  element 

» 

SN(X,D  =  1  /PIVOT  ?  NEXT  IR 
RE  T  URN 
Ap i ot  : 


|  *  * 

!  *  PLOTTING  SUBROUTINE  * 


2570  NLT= 1 

2580  CLEAR  #  IF  NUMPLT=0  Then  GOTO  2780 
2b90  ALPHA  (j.  DISP  "PARAMETERS  OF  FIT"  6  DISP 
2b00  FOR  1=1  TO  NP  9  DISP  AA(I>  $  NEXT  I 
2bl0  IF  NSM  =  2  THEN  GOTO  2640 
2b20  FOR  1=1  TO  NP  9  COOL < I ) =AA( I )  9  NEXT  I 
2b  JO  GOTO  2650 

2b40  FOR  1=1  TO  NP  Li  HEAT ( I )  =AA(  I )  $  NEXT  I 
2650  CHI2=0 

2660  IF  NTEST  <2  THEN  GOTO  2700 
2670  SAR=4.2*(HEAT( 1 > -COOL ( 1  > ) 

2680  IF  NS=2  THEN  SAR=2*SAR*ATIME< JTURN ) +4 . 2* ( HEAT (2>-C00L<2) ) 

2bS0  DISP  FNIVIDSC "SAR  IS" ) , SAR , "Wat ts  per  gram" 

2700  FOR  1=1  TO  N  C*  Y2  =  0  9  FOR  J=1  TO  NS  +  1 
2710  IF  NS+ 1  - J  =  0  AND  BTIME<I>  =  0  THEN  GOTO  2730 
2720  Y2=AA( J)*BTIME(I)  (NS+1 -J)+Y2  9  GOTO  2740 
2/JO  Y2=AA( J)+Y2 
2/40  NEXT  J 

2750  CHI2=<BTEMP(I>- Y2)  2/Y2+CHI2 
2760  NEXT  I 

2770  DISP  "OBSERVED  CHI2=  " .CHI2,”  WITH  ",N*NS+1,"  D.F 

2/80  DISP  "DO  YOU  WANT  A  HARD  COPY  PL01V  Y(£b>  UR  N(Ur 

2790  INPUT  NYSv*  IF  NUMSMO=l  THEN  GOTO  2B9U 

2800  DISP  "TIME  WILL  BE  PLOTTED  FROM  " , XMIN , "TO" , XMAX 

2810  DISP  "TEMPERATURE  WILL  BE  PLOTTED  FROM” , YMIN , "TO" , YMAX 

2820  DISP  "IF  YOU  ARE  SATISFIED  PRESS  RETURN" 

2830  INPUT  NYL IMS 

2840  IF  NYLIMS=""  THEN  GOTO  2890 

2850  DISP  "ENTER  TIME  RANGE  TO  BE  PLOTTED,  TMIN.TMAX" 

2«bU  INPUT  XMIN, XMAX 

2b70  DISP  "ENTER  TEMPERATURE  RANGE  TO  BE  PLOTTED" 

2880  INPUT  YM1N.YMAX 

2890  IF  NTEST>=  1  THEN  GOTO  3050 

2900  CLEAR 

2910  DISP  "CALCULATING  PLOT  SCALING  PARAMETERS  "uFN IV IDS < "PLEASE  WAIT 
2b..u  AAMAX=  In  i  <  lG  I  (  XMAX  >  ) 

2930  A YMAX= I NT  <  LGT  <  YMAX ) ) 

2940  Xr  =10  AXMAX  6  YF  =  1 0  AYMAX 

cj  0  XhWR4  =  VAL.5  (AXMAX)  9  YPWRS=VALS  (AYMAX ) 

29b0  1  XMIN=0  6  YM]N=0 


) 


2970  IF  XMAX/XF>=  1  THEN  DELX=.2 
2980  IF  XMAX/XF>=  3  THEN  DELX-.5 
2990  IF  XMAX/XF>=  5  THEN  DELX=1 
3000  IF  YMAX/YF>=  1  THEN  DELY= .2 
3010  IF  YMAX/YF>=  3  THEN  DELY= . 5 
3020  IF  Yf1AX/YF>=  5  THEN  DEL Y= 1 
3030  XNUN=CEIL  < (XMAX-XMIN)/(XF-DELX) > 

3040  YNUM=CEIL  ( (  YMAX- YMIN) / (  YF *DELY  ) ) 

3050  IF  NY5="N"  THEN  GOTO  3090 
3060  PLOTTER  IS  505 

3070  GRAPH  $  IF  NUMPLT=0  THEN  GCLEAR  §  DEG 
3080  GOTO  3130 
3u90  PLOTTER  IS  1 
3100  GRAPH 

3110  ♦  IF  NUMPLT#0  AND  NTEST>0  THEN  GOTO  3110 

3120  IF  NUMPLT=0  THEN  GCLEAR 

3130  LOCATE  30, RATIO  *100-10,20,95 

3140  !  IF  NUMPLT*0  AND  NTEST>0  THEN  GOTO  3110 

3150  SCALE  XMIN/XF ,XMIN/XF+XNUM*DELX ,  YMIN/YF ,  YMIN/YF+YNUM*DELY 

JlbU  iF  NILST>0  AND  NUMPLl#U  THEN  bU  I U  33bU 

3170  FXD  1 

3180  LGRID  DELX , DELY ,XMIN/XF , YMIN/YF 
3190  AXS=(XMAX-XMIN)/(XF*10> 

3300  AYS=(YMAX-YM1N)/(YF*10) 

3210  MODE  XMIN/XF  +4 , 5*AXS , YMIN/YF -2* AYS  6  LORG  6 
3220  LABEL  USING  "K"  ;  "TIME  (SECS)  x  lE"iXPWR$ 
b23U  MUVt  XMIN/XF -2* AXS, YMIN/YF +4. 5*A)S 
3240  LORG  4  &  DEG  «  LDIR  90 

3250  LABEL  USING  "K"  ;  "TEMP .  (DEGS  C>  x  1E"4YPWRS> 

3260  LDIR  0 

3270  IF  NPRGBE=2  THEN  GOTO  3340 
3280  LINE  TYPE  1 

3290  MOVE  ATIriE< JMIN>/XF , ATEMPC JNIN)/ YF 
33UU  FOR  I-JMIN  TO  JMAX 
3310  DRAW  ATIME(I)/XF,ATEMP(I>/YF 
3320  NEXT  I 

3330  IF  NPROBE= 1  THEN  GOTO  3370 
3340  LINE  TYPE  5 

3350  MOVE  BTIME(KNIN)/XF , BTEMP(KMIN) /YF 

3360  FOR  I-KPIIN  TO  KMAX  e  DRAW  BTIME(  I ) /XF  ,  BTEMP(  I )  /  YF  &  NEXT  I 
3370  IF  NYSLS="N"  OR  NYSL$="NO"  THEN  GOTO  3430 
3380  LINE  TYPE  4 

3390  MOVE  BTIME(801 )/XF,BTEMP(801 )/YF 

3400  FOR  1=80 1  TO  1000 

3410  DRAW  BTIME(I)/XF,BTEMP(I)/YF 

3420  NEXT  I 

3430  PEN  UP 

3440  MOVE  XMIN/XF, YMIN/YF 

3450  IF  SAR=0  THEN  GOTO  3490 

3460  MOVE  XMIN/XF+4.5*AXS, YMAX/YF+ . 1 *AYS 

3470  LORG  4 

3480  LABEL  "SAR” , SAR , "Wat ts  per  gram" 

3490  NTEST-NTEST+1 
3500  PAUSE 
3510  KEY  LABEL 
3520  RETURN 
3530  CHAIN  "DUMP" 

3540  END 
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INTRODUCTION 

Dielectric  heating  is  a  recognized  mechanism  for  the  induction  of  radiofrequency  (RF); 
bioeffects  considered  to  be  thermal  in  origin  -  that  is  attributable  to  temperature  increases. 
However,  apart  from  RF-heat  cell  killing  (e.g.  Sapareto  et  al.,1982,  Chang  et  al., 1987)  very 
little  information  is  available  on  quantitative  relationships  between  RF  bioeffects  at  the 
cellular  and  molecular  level  and  temperature  profiles  over  time  or  thermal  dosage  (TD).  The 
determination  of  RF  heating/cooling  curves  can  be  used  to  compute  the  specific  absorption 
rate  (SAR)  and  the  specific  absorption  (SA)  in  an  exposed  sample  (Stuchly  and  Stuchly, 
1986).  Biological  variables  and  exposure  conditions  can  be  controlled  in  experiments  with  in' 
vitro  systems  to  a  degree  not  achievable  in  vivo.  An  attractive  model  to  study  RF  bioeffects 
is  the  transformation  of  lymphocytes  in  vifro,  provided  the  biological  variables  can  be  related 
to  dosimetric  quantities  that  characterize  exposure  (Czerski,1975;  Budd  and  Czerski,1985). 
To  accomplish  this,  we  designed  an  exposure  system  with  provision  for  real-time  temperature 
monitoring  with  RF-field  non-perturbing  temperature  probes.  The  exposure  system,  which 
will  be  described  in  detail  later,  has  multiple  sample  chambers.  To  allow  on-line  thermometry 
and  dosimetry,  one  of  these  chambers  is  used  as  a  site  for  a  non-perturbing  temperature  probe. 
From  the  temperature  (T)/time(t)  history  of  the  sample  chamber,  an  exposure  dosage  can  be 
determined  and  described  in  terms  of  TD,  SAR,  SA,  or  electric  field  strength.  The  system  can 
be  used  for  studies  of  RF-bioeffects  in  any  tissue  culture  cell  line  or  other  in  vitro  biological ! 
sample.  The  monitoring  component  can  be  applied  to  the  study  of  temperature-dependent 
effects,  irrespective  of  the  modality  used  for  heating. 

EXPOSURE  SYSTEM 

A  schematic  diagram  of  the  system  for  microwave  exposure  of  samples  is  shown  in  Fig.  1. 


The  output  from  a  C W  Hewlett-Packard  Model  8616A  oscillator  feeds  a  Hughes  Model  1 177H 
TWT  amplifier.  Pulsed  or  amplitude-modulated  waveforms  are  obtained  through  the  use  of 
a  Hewlett-Packard  Model  8403A  p-i-n  modulator  and  driver.  For  high  peak  power,  but  low 
duty-cycle,  pulsed  exposures  the  signal  can  be  further  amplified  with  a  IkW  amplifier  (MCL, 
La  Grange,  Ill.)  The  amplified  signal  enters  a  shorted  section  of  S-band  rectangular  waveguide 
300  mm  long  through  a  matched  coaxial  feed.  The  waveguide  sample  holder  receives  its 
power  through  an  isolator  (not  shown  specifically  in  Fig.  l)  and  a  dual- directional  coupler 
that  allows  forward  and  reflected  powers  to  be  measured  with  Hewlett-Packard  Model  432A 
power  meters,  and  waveforms  to  be  monitored.  A  sample  holder  for  biological  specimens 
is  supported  centrally  in  the  waveguide  in  a  block  of  low-density  polystyrene  that  fills  the 
cross-section  of  the  guide.  An  identical  waveguide  section  with  an  identical  plastic  sample 
holder  is  used  for  sham  exposures.  During  microwave  exposures,  both  waveguide  assemblies 
were  housed  in  a  CO2  tissue  culture  incubator  thermostated  at  37. CPC.  Various  sample 
holders  were  tested;  a  4-chamber  plastic  tissue  culture  slide  (Miles  Scientific  #  4804)  was 
found  satisfactory  in  our  application.  Each  chamber  holds  1  ml.  A  miniature  thermistor 
temperature  probe  enters  each  waveguide  section  through  a  small  hole  in  the  top  shorting 
plate.  The  hole  is  drilled  near  the  guide  wall,  at  the  center  of  the  shorter  dimension  -  this  is 
a  low  electric  field  point  and  the  hole  produces  minimal  disturbance.  The  active  end  of  the 
thermistor  probe  enters  the  sample  in  one  of  the  4  chambers  of  the  tissue  culture  dish.  We 
could  move  the  probe  from  chamber  to  chamber  to  check  exposure  uniformity.  Two  types  of 
non-field  perturbing  probe  have  been  used  in  this  way:  A  Narda  Model  801  lB  non-perturbing 
double  temperature  probe  or  two  Vitek  Model  101  probes.  The  temperatures  in  both  the 
exposed  and  sham-exposed  samples  are  recorded  continuously  during  an  experimental  run. 
The  temperature  probes  are  connected  to  a  Hewlett-Packard  relay  activator,  Model  59306A, 
which  is  itself  connected  to  a  Keithley  Model  192  DVM.  Both  the  relay  activator  and  DVM 
are  under  the  control  of  a  Hewlett-Packard  Model  86  desk-top  computer  through  the  HPIB 
(IEEE-488  bus) (Fig. 2). 


Fig. 2.  Schematic  diagram  of  multiple-chamber  sample  holder  with  non¬ 
perturbing  temperature  measurement  and  automated  dosimetry  system. 


Under  computer  control,  the  temperatures  of  exposed  and  sham-exposed  samples  are 
recorded  sequentially  and  stored  in  memory.  Since  the  temperature  of  sham-exposed  samples 
is  essentially  constant,  it  is  monitored  less  frequently  than  the  temperature  of  the  exposed 
samples.  A  typical  experimental  protocol  would  involve  10  temperature  readings  of  the  ex¬ 
posed  sample  for  every  one  reading  of  the  sham-exposed.  For  conventional  thermal  exposures, 
the  temperature  can  be  elevated  by  increasing  the  temperature  of  the  incubator. 

TEMPERATURE/TIME  ANALYSIS 

During  a  typical  experiment,  the  temperature  is  recorded  at  regular  intervals,  with  a 
minimum  measurement  interval  <  la,  before  the  beginning  of,  during,  and  after  the  exposure. 
The  “on”  and  “off”  times  of  exposure  are  recorded  on  the  computer  by  the  operator  using  a 
“soft-key”  interrupt  capability.  At  the  conclusion  of  a  run,  the  T(t)  behavior  is  analyzed  to 
determine  SAR. 

If  the  SAR  is  S  (Wkg~l)  and  the  specific  heat  of  the  sample  is  C(Jkg_1  K~*),  then  the 
SAR  can  be  determined  from  the  change  in  heating  rate  when  microwave  power  is  applied. 


In  practice,  this  analysis  is  performed  automatically.  At  the  end  of  a  run,  the  point  of 
inflection,  or  turning  point,  in  T(t)  is  found  numerically.  The  temperature  profile  to  the  left 
and  right  of  this  point  is  fitted  by  least-squares  to  a  linear  or  quadratic  function  and  the 
change  in  slope  at  the  beginning  of  the  exposure  period  yields  the  SAR.  Figures  3  and  4  are 
examples  of  such  a  procedure. 


Time  jocund*)  *  100 


Fig. 3.  Temperature/time  history  of  an 
exposed  sample  showing  linear  fits  to  cool¬ 
ing  and  heating  portions  of  curve  for  dosime¬ 
try. 
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Fig. 4.  Temperature/time  history  of  ex¬ 
posed  (upper  curve)  and  sham-exposed 
samples  (lower  curve).  The  dosimetry 
for  the  exposed  sample  has  been  deter¬ 
mined  from  a  quadratic  fit  to  the  heating 
curve  during  microwave  exposure  and  a 
linear  fit  to  the  equilibrated  portion  of 
the  curve  prior  to  exposure. 


Figure  3  shows  the  temperature/time  history  of  a  sample  that  was  cooling  prior  to  the  start 
of  microwave  exposure.  The  linear  fits  to  the  cooling  and  heating  portions  of  the  curve  yield 
the  SAR.  Figure  4.  shows  the  heating  of  a  sampled  that  was  quite  well  equilibrated  before 
exposure.  Because  the  signal/noise  ratio  is  high  the  SAR  was  determined  by  a  quadratic 
fit  to  the  heating  portion  of  the  curve.  The  quality  of  the  fit  near  the  turning  point  can 
always  be  examined  to  determine  whether  the  slope  value  is  realistic.  However, unless  the 
temperature/time  profile  has  very  low  noise,  it  is  generally  better  to  use  the  linear  fit  to  T(t) 
near  the  turning  point  or  a  biased  estimate  of  the  slope  can  result.  Figure  4  also  shows  the 
temperature/time  history  of  the  sham-exposed  sample  to  illustrate  its  temperature  stability. 

RESULTS  AND  DISCUSSION 

The  advantages  of  the  procedure  described  above  are  severalfold.  The  sample  need  not 
be  equilibrated  before  SAR  is  determined.  Repeated  determination  of  SAR  can  easily  and 
quickly  be  made,  which  allows  the  SAR  uniformity  from  one  sample  chamber  to  another 
to  be  determined.  Exact  knowledge  of  the  microwave  power  is  not  required:  reproducible 
exposures  can  be  made  at  known  SAR  by  using  the  measured  forward  power  corresponding 
to  a  given  SAR  measurement,  provided  the  experimental  arrangement  is  not  altered  between 
exposures.  Dosimetry  is  not  affected  by  other  losses  in  the  system.  The  waveguide  exposure 
system  with  four  sample  chambers  was  found  to  give  an  SAR  uniformity  from  one  chamber 
to  the  other  within  10%.  Thus,  experimental  samples  can  be  exposed  and  examined  in 
triplicate,  the  fourth  chamber  being  used  for  the  insertion  of  the  temperature  probe. 

The  system  is  biocompatible  and,  depending  on  cell  line,  cell  density,  and  medium,  allows 
continuous  exposures  of  several  days  duration.  The  biocompatibility  of  this  arrangement  was 
proven  by  studying  the  growth  of  human  lymphocytes  under  various  conditions  within  the 
waveguide.  As  well  as  providing  biocompatible  exposure  conditions,  the  system  is  flexible: 
it  allows  the  exposure  of  tissue  culture  cells  growing  in  suspension  or  in  monolayers,  and  can 
easily  be  used  with  cell-free  samples.  The  exposure  system  is  relatively  simple,  and  can  be 
assembled  from  off-the-shelf  components.  Its  capabilities  can  be  expanded  by  introducing 
computer  control  of  RF  power  input  based  on  feedback  from  temperature  measurements. 
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ABSTRACT 


A  waveguide  exposure  system  with  automated  sample  temperature  measurement  is  de¬ 
scribed,  This  system  provides  on-line  determination  of  the  thermal  dose  received  by  biological 
samples  in  vitro.  It  allows  automated  computation  of  the  specific  absorption  rate  determined 
from  heating/cooling  curves,  uses  non-perturbing  thermometry,  is  biocompatible,  and  can  be 
used  for  measurements  of  both  microwave  and  conventional  heating. 


INTRODUCTION 


Dielectric  heating  is  a  recognized  mechanism  for  the  induction  of  radiofrequency  (RF) 
bioeffects  considered  to  be  thermal  in  origin  -  that  is  attributable  to  temperature  increases. 
However,  apart  from  RF-heat  cell  killing  (e.g.  Sapareto  et  ah, 1982,  Chang  et  al.,1987)  very  lit¬ 
tle  information  is  available  on  quantitative  relationships  between  RF  bioeffects  at  the  cellular 
and  molecular  level  and  temperature  profiles  over  time  or  thermal  dosage  (TD).  The  deter¬ 
mination  of  RF  heating/cooling  curves  can  be  used  to  compute  the  specific  absorption  rate 
(SAR)  and  the  specific  absorption  (SA)  in  an  exposed  sample  (Stuchly  and  Stuchly,  1986). 
Biological  variables  and  exposure  conditions  can  be  controlled  in  experiments  with  in  vitro 
systems  to  a  degree  not  achievable  in  vivo.  An  attractive  model  to  study  RF  bioeffects  is 
the  transformation  of  lymphocytes  in  vitro,  provided  the  biological  variables  can  be  related 
to  dosimetric  quantities  that  characterize  exposure  (Czerski,1975;  Budd  and  Czerski,1985) . 
To  accomplish  this,  we  designed  an  exposure  system  with  provision  for  real-time  temperature 
monitoring  w'ith  RF-field  non-perturbing  temperature  probes.  The  exposure  system,  which 
will  be  described  in  detail  later,  has  multiple  sample  chambers.  To  allow  on-line  thermometry 
and  dosimetry,  one  of  these  chambers  is  used  as  a  site  for  a  non-perturbing  temperature  probe. 
From  the  temperature  (T)/time(t)  h:  tory  of  the  sample  chamber,  an  exposure  dosage  can  be 
determined  and  described  in  terms  of  TD,  SAR,  SA,  or  electric  field  strength.  The  system  can 
be  used  for  studies  of  RF-bioeffects  in  any  tissue  culture  cell  line  or  other  in  vitro  biological 
sample.  The  monitoring  component  can  be  applied  to  the  study  of  temperature-dependent 
effects,  irrespecthe  of  the  modality  used  for  heating. 

EXPOSURE  SYSTEM 

A  schematic  diagram  of  the  system  for  microwave  exposure  of  samples  is  shown  in  Fig. 
1  The  output  from  a  CW  Hewlett-Packard  Model  8616A  oscillator  feeds  a  Hughes  Model 
1177H  TWT  amplifier.  Pulsed  or  amplitude-modulated  waveforms  are  obtained  through  the 
use  of  a  Hewlett-Packard  Model  8403A  p-i-n  modulator  and  driver.  For  high  peak  power,  but 
low  duty-cycle,  pulsed  exposures  the  signal  can  be  further  amplified  with  a  lkW  amplifier 
(MCL,  La  Grange,  Ill.)  The  amplified  signal  enters  a  shorted  section  of  S-band  rectangular 
waveguide  300  mm  long  through  a  matched  coaxial  feed.  The  waveguide  sample  holder  receives 
its  power  through  an  isolator  (not  shown  specifically  in  Fig,  l)  and  a  dual-direct, onal  coupler 
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that  allows  forward  and  reflected  powers  to  be  measured  with  Hewlett-Packard  Model  432A 
power  meters,  and  waveforms  to  be  monitored.  A  sample  holder  for  biological  specimens  is 
supported  centrally  in  the  waveguide  in  a  block  of  low-density  polystyrene  that  fills  the  cross- 
section  of  the  guide.  An  identical  waveguide  section  writh  an  identical  plastic  sample  holder 
is  used  for  sham  exposures.  During  microwave  exposures,  both  waveguide  assemblies  were 
housed  in  a  CO2  tissue  culture  incubator  thermostated  at  37. CPC.  Various  sample  holders 
were  tested:  a  4-chamber  plastic  tissue  culture  slide  (Miles  Scientific  #  4804)  was  found 
satisfactory  in  our  application  (Fig. 2).  Each  chamber  holds  1  ml.  A  miniature  thermistor 
temperature  probe  enters  each  waveguide  section  through  a  small  hole  in  the  top  shorting 
plate.  The  hole  is  drilled  near  the  guide  wall,  at  the  center  of  the  shorter  dimension  -  this  is 
a  low  electric  field  point  and  the  hole  produces  minimal  disturbance.  The  active  end  of  the 
thermistor  probe  enters  the  sample  in  one  of  the  4  chambers  of  the  tissue  culture  dish.  We 
could  move  the  probe  from  chamber  to  chamber  to  check  exposure  uniformity.  Two  types  of 
non-field  perturbing  probe  have  been  used  in  this  way:  A  Narda  Model  801  IB  non-perturbing 
double  temperature  probe  or  two  Vitek  Model  101  probes.  The  temperatures  in  both  the 
exposed  and  sham-exposed  samples  are  recorded  continuously  during  an  experimental  run. 
The  temperature  probes  are  connected  to  a  Hewlett-Packard  relay  activator,  Model  59306A, 
which  is  itself  connected  to  a  Keithley  Model  192  DVM.  Both  the  relay  activator  and  DVM 
are  under  the  control  of  a  Hewlett-Packard  Model  86  desk-top  computer  through  the  HPIB 
(IEEE-488  bus).  Under  computer  control,  the  temperatures  of  exposed  and  sham-exposed 
samples  are  recorded  sequentially  and  stored  in  memory.  Since  the  temperature  of  sham- 
exposed  samples  is  essentially  constant,  it  is  monitored  less  frequently  than  the  temperature 
of  the  exposed  samples.  A  typical  experimental  protocol  would  involve  10  temperature  readings 
of  the  exposed  sample  for  every  one  reading  of  the  sham-exposed.  For  conventional  thermal 
exposures,  the  temperature  can  be  elevated  by  increasing  the  temperature  of  the  incubator. 

TEMPERATURE/TIME  ANALYSIS 

During  a  typical  experiment,  the  temperature  is  recorded  at  regular  intervals,  with  a 
minimum  measurement  inter  al  <  Is,  before  the  beginning  of,  during,  and  after  the  exposure. 
The  “on*  and  “ofT  times  of  exposure  are  recorded  on  the  computer  by  the  operator  using  a 
“soft-key"  interrupt  capability.  At  the  conclusion  of  a  run,  the  T(t)  behavior  is  analyzed  to 
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determine  SAR. 


If  the  SAR  is  5  (Wkg-1)  and  the  specific  heat  of  the  sample  is  C(Jkg“ 1  K  ~ 1 ),  then  the 
rate  of  heating  during  microwave  exposure,  for  small  temperature  variations  over  which  C  can 
be  assumed  to  be  constant,  is 


If,  prior  to  deliberate  heating,  a  sample  is  below  its  equilibrium  temperature,  Tcq,  with  its 
surroundings,  then  for  small  temperature  differences  from  equilibrium,  the  rate  of  temperature 
change  is  of  the  form 

(§)+  =  *(r„-T).  (2) 

Similarly,  if  after  an  exposure  is  terminated  a  sample  cools  the  rate  of  temperature  change 
will  be 

(?)  =  ~a(T‘*  -  T)-  (3) 

Eqs.  (2)  and  (3)  are  equivalent  statements  of  Newton's  law  of  cooling.  In  general,  when  mi* 
crowave  power  is  applied  to  a  sample  that  was  not  in  thermal  equilibrium  with  its  surroundings 
at  the  start  of  exposure: 

dT  S 

~di  =  c  +  “(T"  -  T)  (<) 

The  solution  to  Eq.  (4)  can  be  written  in  the  form 

I'-T”=(sl  +  r--:r‘)<1  -•■*’)  (5) 

where  Ta  is  the  control  temperature  at  t  =  0.  Either  Eq.  (4)  or  Eq.  (5)  can  be  used  to 
determine  the  SAR.  From  Eq.  (4), 

so  that  determination  of  the  heating  rate  before  and  during  application  of  microwave  power, 
or  during  and  after  the  application  of  microwave  power,  determines  the  SAR. 
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If  the  sample  is  in  thermal  equilibrium  before  exposure  begins,  Eq.  (5)  reduces  to  the 
simpler  form 

T~To  =  ^{1~ea<)  (7) 

which  can  be  used  to  find  S  from  T(t).  In  general,  the  approach  using  Eq.  (6)  is  preferable, 
since  for  times  near  t  —  0,  the  exponential  can  be  expanded  in  quadratic  form. 


In  practice,  this  analysis  is  performed  automatically.  At  the  end  of  a  run,  the  point  of 
inflection,  or  turning  point,  in  T(t)  is  found  numerically.  The  temperature  profile  to  the  left 
and  right  of  this  point  is  fitted  by  least-squares  to  a  linear  or  quadratic  function  and  the  change 
in  slope  at  the  beginning  of  the  exposure  period  yields  the  SAR.  Figures  3  and  4  are  examples 
of  such  a  procedure.  Figure  3  shows  the  temperature/time  history  of  a  sample  that  was  cooling 
prior  to  the  start  of  microwave  exposure.  The  linear  fits  to  the  cooling  and  heating  portions 
of  the  curve  yield  the  SAR.  Figure  4.  shows  the  heating  of  a  sample  that  was  quite  well 
equilibrated  before  exposure.  Because  the  signal/noise  ratio  is  high  the  SAR  was  determined 
by  a  quadratic  fit  to  the  heating  portion  of  the  curve.  The  quality  of  the  fit  near  the  turning 
point  can  always  be  examined  to  determine  whether  the  slope  value  is  realistic.  However, unless 
the  temperature/time  profile  has  very  low  noise,  it  is  generally  better  to  use  the  linear  fit  to 
T(t)  near  the  turning  point  or  a  biased  estimate  of  the  slope  can  result.  Figure  4  also  shows  the 
temperature/time  history  of  the  sham-exposed  sample  to  illustrate  its  temperature  stability. 

RESULTS  AND  DISCUSSION 


The  advantages  of  the  procedure  described  above  are  severalfold.  The  sample  need  not 
be  equilibrated  before  SAR  is  determined.  Repeated  determination  of  SAR  can  easily  and 
quickly  be  made,  which  allows  the  SAR  uniformity  from  one  sample  chamber  to  another  to  be 


determined.  Exact  knowledge  of  the  microwave  power  is  not  required:  reproducible  exposures 
can  be  made  at  known  SAR  by  using  the  measured  forward  power  corresponding  to  a  given 
SAR  measurement,  provided  the  experimental  arrangement  is  not  altered  between  exposures. 
Dosimetry  is  not  affected  by  other  losses  in  the  system.  The  waveguide  exposure  system  with 
four  sample  chambers  was  found  to  give  an  SAR  uniformity  from  one  chamber  to  the  other 


within  10%.  Thus,  experimental  samples  can  be  exposed  and  examined  ,n  triplicate,  the  fourth 
chamber  being  used  for  the  insertion  of  the  temperature  probe. 
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The  system  is  biocompatible  and,  depending  on  cell  line,  cell  density,  and  medium,  allows 
continuous  exposures  of  several  days  duration.  The  biocompatibility  of  this  arrangement  was 
proven  by  studying  the  growth  of  human  lymphocytes  under  various  conditions  within  the 
waveguide.  106  cells/ml  in  chromosome  medium  1 A  (Gibco)  with  or  without  phytohemagglutin 
were  incubated  for  72  or  120  hours  in  the  exposure  system  wuth  no  applied  RF  power.  Cell 
morphology,  lymphoblastoid  transformation,  mitotic  index,  and  cell  viability  tested  by  neutral 
red,  Janus  green  and  trypan  blue  stains  did  not  differ  from  those  in  lympocyte  cultures  from 
the  same  donor  cultured  in  a  conventional  CO2  incubator  according  to  the  protocol  provided 
by  Gibco  with  their  media.  As  well  as  providing  biocompatible  exposure  conditions,  the  system 
is  flexible:  it  allows  the  exposure  of  tissue  culture  cells  growing  in  suspension  or  in  1  molayers, 
and  can  easily  be  used  with  cell-free  samples. 

The  exposure  system  is  relatively  simple,  and  can  be  assembled  from  off-the-shelf  compo¬ 
nents.  Its  capabilities  can  be  expanded  by  introducing  computer  control  of  RF  power  input 
based  on  feedback  from  temperature  measurements,  in  a  similar  way  to  the  system  described 
by  Chang  et  al.(l987).  This  system  used  a  stripline  operating  at  915  MHz,  and  allowed  uniform 
exposure  of  10  ml  tissue  culture  samples.  The  temperature/time  profile  of  heating  could  be 
pre-programmed,  and  biocompatibilty  allowed  exposures  of  several  days  duration.  Guy  (1977) 
designed  a  system  for  exposure  of  5  ml  samples  in  the  range  from  dc  to  1  GHz.  Field  strength 
dosimetry  was  obtained  from  on-line  measurement  of  feed-line  impedance.  The  temperature 
of  the  sample  could  be  maintained  at  a  preset  level  using  a  circulating  liquid  heat  exchanger. 
No  data  on  biocompatability  of  this  system  were  provided,  however  the  information  given  sug¬ 
gests  suitability  only  for  short-term  exposures.  Lin(l976)  proposed  the  use  of  micropipettes 
immersed  in  a  fluid-filled  waveguide  irradiation  hamber  for  in  vitro  studies  under  controlled 
temperature  and  dosimetric  conditions.  This  system  was  also  only  suitable  for  short-term 
exposures.  Allis  et  al.  (1975)  installed  a  waveguide  exposure  system  in  a  Cary  15  spectropho¬ 
tometer  allowing  simultaneous  exposure  and  spectrophotometric  biochemical  determinations. 
Our  system  and  the  ones  described  in  the  papers  quoted  here  provide  a  wide  range  of  pos¬ 
sibilities  for  in  vitro  studies  under  controlled  conditions,  including  precise  RF  dosimetry.  By 
using  one  of  these  systems,  or  combining  desired  features  of  some  of  them,  most  requirements 
for  the  study  of  RF  bioeffects  at  the  cellular  and/or  molecular  level  can  be  met. 
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FIGURE  CAPTIONS 

1.  Schematic  diagram  of  essential  components  of  microwave  exposure  system. 

2.  Schematic  diagram  of  multiple-chamber  sample  holder  with  non-perturbing  temperature 
measurement  and  automated  dosimetry  system. 

3.  Temperature/time  history  of  an  exposed  sample  showing  linear  fits  to  cooling  and  heating 
portions  of  curve  for  dosimetry. 

4.  Temperature/time  history  of  exposed  (upper  curve)  and  sham-exposed  samples  (lower 
curve).  The  dosimetry  for  the  exposed  sample  has  been  determined  from  a  quadratic  fit  to 
the  heating  curve  during  microwave  exposure  and  a  linear  fit  to  the  equilibrated  portion  of 
the  curve  prior  to  exposure. 
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